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THE  APPLICATION  OF  OCCULTATIONS  TO  GEODESY 


PART  1:  GENERAL  THEORY 


Introduction 


The  sun  and  the  moon  have  apparent  daily  motions  across 
the  sky,  rising  in  the  east  and  setting  in  the  west.  These 
motions  are  apparent  in  that  they  are  presumably  not  motions 
of  the  sun  and  moon  proper,  but  are  the  result  of  the  earth's 
rotation.  Hence  these  motions  are  shared  by  all  the  rest  of 
the  heavenly  bodies  —  the  stars  and  the  planets.  That  the 
sun  and  moon  have  other  motions  which  are  not  related  to  the 
earth's  rotation  is  shown  by  the  different  positions  of  these 
two  from  day  to  day  with  respect  to  the  stars.  This  is 
especially  obvious  in  the  case  of  the  moon,  which  can  be  seen 
to  shift  its  position  eastward  with  respect  to  the  stars  by  a 
distance  equal  to  its  own  diameter  every  hour.  Gince  the 
moon  has  an  angular  diameter  of  about  half  a  decree,  in  one 
day  of  2b  hours  the  moon  will  have  moved  12°  eastward  with 
respect  to  the  background  stars,  and  hence  will  be  observed  to 
rise  the  time  equivalent  of  12°,  or  b8  minutes,  later  every 
day.  A  complete  circuit  of  the  star  field  takes  about  27.3 
days. 


In  its  eastward  motion  among  the  stars,  it  is  to  be 
expected  that  the  mow  will  pass  in  front  of  and  hide  such 
luminous  bodies  as  may  be  along  its  path.  Except  where  the 
body  is  our  own  sun,  in  which  case  an  eclipse  is  said  to 
occur,  the  passage  of  the  moon  across  the  line  of  sight 
between  a  star  or  planet  and  an  observer  is  called  an 
occultation  of  the  star  or  planet  by  the  moon.  An  eclipse 
of  one  sort  or  another  seldom  occurs  more  than  twice  a  year; 
an  occultation,  on  the  other  hand,  can  be  observed  almost 
every  night  with  a  pair  of  binoculars. 

Both  eclipses  and  occultations  have  been  observed  in 
order  to  find  the  geographic  location  of  the  observer  (ref¬ 
erences  12,  13,  lb).  That  such  use  can  be  made  of  eclipses 
and  occultations  stems  basically  from  three  circumstances: 

(1)  The  moon  is  so  close  to  the  earth  that  a 
shift  in  the  observer's  position  produces  an  appreciable 
change  in  the  apparent  position  of  the  moon  with  respect 
to  the  surrounding  stars  (or  the  sun).  At  the  sub -lunar 
point,  a  shift  of  1000  meters  in  the  observer's  position 
produces  a  shift  of  about  one-half  second  of  arc  in  the 
position  of  the  moon. 
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(2)  The  directions  of  the  brighter  stars  occulted 
by  the  moon  are  well  known,  so  that  a  bright  star  occultation 
gives  an  accurate  space  direction  for  that  portion  of  the 
edge  of  the  moon  (at  the  instant  of  occultation)  at  which  the 
occultation  occurs. 

(3)  The  position  of  the  moon  with  respect  to  the 
earth  is  known  from  the  exhaustive  studies  of  E.  W.  Brown 
(reference  1)  and  others,  as  a  result  of  many  years  of  precise 
observations  at  observatories.  At  the  instant  of  occultation 
there  is  established  in  space  a  fixed  line  on  which  the  observer 
is  known  to  lie.  Since  the  occultation  takes  place  at  the  edge 
of  the  moon  and  not  at  the  center,  to  fix  the  line  would  require 
not  only  the  coordinates  of  the  center  of  the  moon,  but  also  the 
direction  and  distance  from  the  center  of  the  moon  to  that  point 
on  the  edge  at  which  the  occultation  was  observed  to  take  place 
(Figure  1). 

In  practice,  the  direction  (from  the  moon's  center  to  its 
edge)  at  which  the  occultation  occurs  cannot  be  measured 
accurately,  and  no  attempt  is  made  to  measure  it,  so  that  there 
is  actually  a  set  of  parallel  lines  determined,  all  pointing 
toward  the  same  star  and  on  any  one  of  which  the  observer  could 
lie.  The  intersection  of  this  set  of  lines  with  the  spheroid 
(taken  to  be  the  earth's  surface)  is  an  arc,  a  short  segment  of 
which  is  called  a  "line  of  position* (LOP).  The  direction  from 
the  center  is  always  known  well  enough  that  on  a  map  the  LOP  is 
practically  a  straight  line.  Two  occultations  give  two  lines 
of  position,  and  the  observation  site  lies  at  the  intersection 
of  these  two  lines. 

Another  deviation  of  practice  from  theory  occurs  because 
the  lunar  profile  is  not  a  perfect  circle,  but  is  highly  irregular; 
that  point  on  the  edge  of  the  moon  at  which  the  occultation  occurs 
is  at  a  different  distance  from  the  lunar  center  for  each  occulta- 
tion  and  this  distance  is  not  accurately  known*  Hence,  every 
occultation  occurs  at  a  different  distance  from  the  optical  center 
of  the  moon;  every  additional  occultation  introduces  another  and 
different  radius  of  the  moon,  and  the  ensuing  set  of  equations 
always  contains  more  unknowns  than  there  are  knowns.  To  remove 
this  difficulty  each  occultation  can  be  observed  by  two  different 
observers;  if  the  occultation  occurs  at  the  same  point  on  the  limb 
for  b  th,  one  of  these  observations  may  be  used  to  determine  the 
radius  of  the  moon  for  that  occultation.  When  this  is  done,  the 
situation  is  like  that  in  the  preceding  paragraph  where  only  the 
position  angle  (direction  from  the  moon's  center)  in  addition  to 
the  observer's  position,  was  unknown,  and  the  problem  is  again 
solvable. 
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FIGURE  1 


The  star  at  effectively  infinite  distance,  is  the  center  for  projection  of -the 
profile  of  an  irregularly-shaped  body,  the  Moon,  onto  an  irregularly-shaped 
surface,  the  Earth.  The  Moon's  surface  is  approximated  by  a  sphere,  its 
profile  by’a  circle  and  the  Earth's  surface  by  an  ellipsoid  of  revolution. 


I .  Mathematical  Formulation  of  tha  Problem 


Mathematically,  tha  a lap last  approach  to  tha  problem 
would  ba  to  go  directly  to  the  solution  by  determining  the 
intersection  of  tha  spheroid  surface  with  tha  set  of  lines 
tangent  to  tha  moon's  adgs  and  extending  in  the  direction  of  the 
occulted  star.  Such  an  approach  would,  however,  be  very  diffi¬ 
cult  to  handle,  since  it  would  entail  the  solution  of  a  complex 
system  of  simultaneous  fourth-degree  aquations.  Hence  a  less 
direct  t;oute  to  the  solution  was  adopted  long  ago  by  astronomers 
and  has  been  retained  ever  since.  This  Indirect  route  permits 
most  of  the  computations  to  be  carried  out  as  if  the  occultation 
occurred  in  a  plane,  and  then  brings  the  results  from  the  plane 
up  to  the  surface  of  the  spheroid  ss  a  last  step. 

As  was  mentioned  in  the  Introduction,  the  two  major  unknowns, 
other  than  the  geodetic  ones,  occurring  in  the  occultation  problem 
are  the  position  angle  to  the  occultation  and  the  distance  from 
the  center  of  the  moon  to  the  point  on  the  edge  at  which  the 
occultation  occurs.  After  the  basic  occultation  equations  have 
been  set  up,  as  in  section  B  below,  two  distinct  ways  of  removing 
these  two  unknowns  are  possible.  One,  which  is  called  the  Single- 
Site  Method,  requires  that  occultatlons  be  observed  from  one  site 
over  as  great  a  range  of  position  angles  as  possible.  When  the 
geodetic  unknowns  and  as  many  as  possible  of  the  systematic 
astronomical  unknowns  are  solved  for  by  the  method  of  least 
squares,  and  the  variations  in  the  radius  of  the  moon  are  ignored, 
it  is  expected  that  variations  in  the  moon's  radius  will  act  in 
the  solution  as  if  they  were  randomly  distributed  with  respect 
to  the  various  observations , and  hence  will  not  appear  in  the 
values  found  for  the  geodetic  unknowns.  This  method  is  discussed 
only  briefly  in  section  C  below,  because  although  it  is  of  con¬ 
siderable  Interest  in  itr  own  right,  as  well  as  closely  related  in 
theory  to  the  Equal -Limb -Line  Method,  it  has  already  been  treated 
at  length  in  a  previous  paper  by  Famelia  A.  Henriksen  (reference  2). 
It  is  included  here  largely  for  the  sake  of  completeness  and  to 
throw  additional  light  on  the  less  obvious  problems  associated 
with  occultation  survey  by  the  Equal -Limb -Line  Method. 

The  Equal -Limb-Line  Method,  as  stated  in  the  'introduction, 
removes  the  uncertainty  in  the  radius  of  the  moon  by  using  a  second 
observer  to  determine  the  radius  for  the  same  occultation  as  that 
being  viewed  by  the  first  observer.  Since  no  reliance  is  placed 
on  statistics  to  produce  results,  the  accuracy  of  this  method  is 
inherently  greater  than  that  of  the  Single-Site  Method.  Section 
D  below  covers  the  Equal-Limb-Line  Method  in  detail. 
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FIGURE  2 

Geodetic  and  topocentric  coordinates  of  the  observer  (x,$,h,  and  E 1  systems) 
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Before  anything  is  done  toward  setting  up  the  equations 
describing  the  occultation  phenomenon,  however,  a  clear  statement 
must  be  made  concerning  the  frames  of  reference  in  which  the 
equations  will  be  used.  This  is  not  a  simple  subject,  since  the 
data  which  come  from  the  observations ,  from  the  epheaaris  of  the 
moon,  and  from  the  star  catalogs  are  referred  to  many  distinct 
frames  of  reference  which  are  related  to  each  other  in  a  complex 
fashion.  Also,  because  the  computation  procedure  is  indirect,  a 
chain  of  intermediate  reference  systems  is  introduced  which  must 
be  distinguished  from  the  primary  reference  systems.  Section  A 
below  will  discuss  simply  the  frames  of  reference  used  in  the 
theory. 

A.  Frames  of  Reference. 

1.  The  observer's  coordinates  (Figure  2}  are  given 
first  with  respect  to  some  standard  spheroid  as: 

a,  h. 

The  sphetoid  is  assumed  to  be  defined  (Boatford,  reference  3) 
so  that  its  minor  axis  is  parallel  to  the  earth's  axis  of 
rotation;  a  major  axis  is  assigned  the  value  A*  0,  and  is 
parallel  to  the  meridian  of  Greenwich , and  the  perpendicular  at 
the  datum  point  7\0,  0O  is  defined  to  coincide  with  the  vertical 
(the  direction  ofgravity)  there  at  a  depth  h^  below  the  surface 
of  the  earth.  Usually,  h  will  be  given  initially  not  with 
respect  to  the  reference  spheroid  but  with  respect  to  mean  sea 
level  —  essentially,  the  co-geold.  The  height  of  the  co-geoid 
above  the  spheroid,  Ah.  must  be  added  to  the  given  hg  to  get 
the  height,  h,  above  the  spheroid. 

Most  computations  are  carried  out  using  the  International 
spheroid  (as  6,378,388  meters,  1/f  a  297)as  the  reference 
surface.  When,  as  often  happens,  the  survey  data  on  the  occultation 
site  are  with  reference  to  some  other  reference  surface,  a  transfor¬ 
mation  between  the  two  spheroids  is  necessary.  The  procedure  used 
is  given  in  Appendix  D. 

2.  A  rectangular  coordinate  system  with  origin  at 

the  center  of  the  spheroid  is  used  in  preference  to  the  geographic 
system  for  computations  (Figure  3);  this  system  is  defined  by  the 
equations : 

X^  ■  (W-h)  cos  A  cos  0 

X2  ■  (N+h)  sin  ^  cos  0  (1) 

x3  -  $(l-e2)  +  Sin  0. 
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N  Is  the  radius  of  curvature  in  the  prime  vertical  and  can  be 
computed  from 

N  -  a  [l-e2  sin2  0  ]  (2) 

a  is  the  semi -major  axis  of  the  spheroid. 


3.  A  topocentric  rectangular  coordinate  system  with 
origin  near  or  at  the  observer  is  used  to  relate  the  occultation 
data  directly  to  the  geodetic  coordinates  of  the  observer.  This 
system  (Figure  2),  designated  the  B*-system.  is  related  to  the 
X*- -system  by  the  equations: 


«* <x‘  - 


sin  X 
cos  ^ 
0 


o 

+  cos  ^ 
-  sin  A 
0 


0 

0 

1 


(3) 


(4) 


1  0  0 

0  +  sin  0  -  cos  0  I  (5) 

L0  +  cos  0  +  sin  0  J 

The  X1  are  the  X* -system  coordinates  of  the  topocentric  system  origin. 


4.  Host  of  the  computations  are  carried  out  in  an  inter¬ 
mediate  reference  system  whose  origin  is  also  at  the  center  of  the 
spheroid  (see  below),  but  whose  orientation  is  with  reference  to  the 
moon's  and  star's  coordinates  rather  than  with  reference  to  earth- 
fixed  coordinates  (references  11.  12,  13).  This  system  will  be  denoted 
by  lower-case  letters  (xl,  x2,  x^).  The  x3-axis  is  parallel  to  the 
line  joining  the  star  and  the  center  of  the  moon  and  is  positive  in  the 
direction  of  the  star;  the  x2-axis  lies  in  the  plane  of  x^  and  the 
spheroid'6  minor-axis  and  is  positive  to  the  North;  the  xl-axis  is 
perpendicular  to  x3  and  is  positive  towards  the  East  (Figure  4).  The 
x-system  is  therefore  related  to  the  Xl-system  by  the  equations: 


[x1]  -  [R  S*1  ]  •  (6) 
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Hera 


M* 


4-  «ln  4-  cos  JLA* 

"  co*  j**  +  #ltt 
0  0 


0 

0 

1 


1 

0 

0 


4*  In  4-  coa 

-cos  4-  aln 


la  the  Greenwich  hour  angle  of  the  ater. 
3  GST  -  flt. 


GST  »  GSTho  4  (tlae  of  obaervatlon) 


(7) 


(8) 


(9) 


If  the  lx7-ayatea  la  referred  to  aome  space-fixed  system 
instead  of  to  an  earth-fixed  system,  the  relation  between  the  two 
can  be  expressed  as  a  rotation: 

coa  0  sin  9  0 

-  sin  9  coa  9  0 

0  0  1 


Here  the  £*X  J  system  has  rotated  through  the  angle  9  since 
zero  time  when  it  made  the  angle  80  with  the  space-fixed  systems, 
and^e,  <f*  are  the  right  ascension  and  declination  of  the  star 
(GST  is  the  Greenwich  hour  angle  of  the  ater  at  the  instant 
considered) . 

In  observing  an  occultation,  the  optical  axis  of  the 
telescope  is  kept  pointed  at  the  star,  and  the  time  at  which 
the  star  disappears  (or  reappears  if  an  emersion  is  being  observed) 
is  recorded.  The  optical  axis  is  therefore  constantly  parallel  to 
to  the  x^-axis,  and  no  distinction  is  made  between  events  occurring 
at  different  points  along  this  axis.  As  far  as  the  observer  is 
concerned,  his  world  of  events  lies  in  the  x1  xz  plane.  Hence  it 
is  customary  to  further  simplify  the  mathematics  by  doing  most  of 
the  computations  in  the  xx  xz  plane,  which  is  called  the 
"fundamental  plane"  after  Bessel,  who  first  introduced  it. 
Coordinates  in  this  plane  era  often  denoted  by  £x,  £z, 
or  | ,  rj . 


fa]  - 
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These  four  systems  are  basic  to  the  mathematics  of  occultatlons. 
There  are  a  few  others,  principally  those  associated  with  the  motion 
of  the  moon  and  with  the  lunar  surface,  but  except  for  the  above  four, 
none  plays  any  great  role  in  the  reduction  of  the  observation  data. 
Considerable  care  must  be  taken,  however,  in  the  use  of  these  four 
systems.  In  the  first  place,  the  observations  give  data  which  are  or 
can  be  referred  to  systems  with  origin  at  the  center  of  the  reference 
spheroid.  The  equations  of  the  moon's  motion,  on  the  other  hand,  are 
set  up  in  a  geocentric  system  --  i.e,,  a  system  with  origin  at  the 
earth's  center  of  gravity;  at  the^present  time,  the  position  of  the 
earth's  center  of  gravity  with  respect  to  any  reference  spheroid  is  . 
known  to  .better  than  ±  100  meters.  Hence  a  certain  amount  of  error 
is  Introduced;  a  more  detailed  study  of  this  point  is  made  in  Part  II, 
In  the  second  place,  if  more  than  one  observing  site  is  used,  and  if 
the  sites  are  on  different  datums,  there  will  be  a  separate  x^-,  X^-, 
E^-system  for  each  site.  The  systems  for  the  separate  sites  will  be 
designated  by  subscripts. 

B.  Basic  Equations 
1.  General: 

There  are  four  physical  bodies  which  enter  into  the 
occultatlon  equations;  these  are  the  star,  the  moon,  the  earth,  and 
the  observer.  They  are  all  tied  together  by  one  parameter  —  the  time. 
The  star's  position  is  given  by  the  vector: 

"r*  -  (  &*.  **). 

where  <**,  ^are  the  star's  apparent  place  coordinates,  but  since 
t^e  distance  r*  is  of  no  Interest  to  occultatlons,  the  unit  vector 
r*(^*.  Sic) ,  which  gives  the  direction  of  the  star,  can  be  used 
Instead.  r+  is  gotten  in  the  usual  manner  from  a  catalog  of  star 
directions;  the  occulted  stars  lie  within  30*  of  the  celestial 
equator.  In  the  above  vectors,  G**,  **  and  r*  are  all  functions  of 
the  time  t,  but  since  the  period  of  time  over  which  a  particular 
occultatlon  can  be  usefully  observed  is  less  than  three  hours 
(usually  less  than  1.5  hours),  r#(t)  for  such  stars  is  practically 
constant. 


The  moon's  position  is  given  by  the  vector: 

where  &m,  &m  and  rB  are  also  time  dependent.  The  elements  of  the 
vector  are  most  easily  gotten  from  one  of  the  published  lunar 
ephemerldes  (reference  6).  They,  of  course,  change  appreciably  within 
a  period  of  three  hours,  and  could  be  written  explicitly  as 

rn  [>(t>’  S*M’ 
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or,  since  ephemerides  give  the  horizontal  equatorial  parallax 
iaatead  of  the  distance,  aa 

r«  j\< t).  ^(t).TTa(t)]. 

In  the  ephemaridea  the  elementa  given  refer  to  a  point  which  ia 
called  the  center  of  maea  of  the  moon,  and  not  to  the  center  of 
the  visible  diak.  Hence  account  must  be  taken,  in  correlating 
obaervatione  with  the  ephemeris  predictions,  of  the  displacement 
of  the  optical  center  from  the  aphemerie  center.  The  equation 
of  a  spheroid  in  a  apace-fixed  ays tern  is 

t  •  S  •  r  -  1  (10) 

where  8  is  a  dyadic  with  elements  which  are  functions  of  time. 

Since  the  occultatlon  occurs  at  the  edge  of  the  moon,  as 
viewed  from  the  earth,  the  equation  of  the  line  on  which  the 
observer  and  star  lie  is: 

**•*» +?«+**  <u> 

where  fm  ia  the  vector  from  the  center  of  the  moon  to  the 
point  on  the  surface  at  which  the  occultatlon  occurs,  and*  is 
a  variable.  The  vector  fL  la  a  complicated  function  of  the 
time  since  the  moon's  motion  ia  compounded  from  many  inde¬ 
pendent  motions  such  as  libratlons,  etc. 

fi  */’.[/.(').  VO,  fmM I  •  (12) 

Here  ^,(t),  bB  (t)-,  P  { t)  are  the  aelenographlc  coordinates  of 
that  point  on  the  moon  at  which  the  line  observer-to-star  af  R 
la  tangent  when  the  occultatlon  occurs.  Our  knowledge  of  the ■ 
shape  of  the  surface  is  insufficient  to  permit  prediction  of/m, 
ba,  and  fm,  and  as  mentioned  before,  certain  assumptions  or 
procedures  must  be  adopted  to  get  around  this  difficulty.  Taking 
it  for  granted  that  this  has  been  done,  the  solution  is  gotten 
by  solving  simultaneously  the  equations: 

/R  -?*  +  ?«+  r*  (13) 

■> 

*  1  -  R  *  S  •  R  (14) 
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Because  of  the  Indeterminacy  in  rm,  Che  first  equation  will 
be  not  that  of  a  single  line,  but  will  define  a  set  of  gener¬ 
ating  elements  of  a  cylinder;  the  common  solution  of  these  two 
equations  will  give  a  curve  instead  of  a  point.  If  n  occupations 
are  observed,  the  vectors  in  these  equations  become  3  x  n  matrices, 
and  we  have 


*[*]  -  [r»]  +  [/V]  +  [r„] 

2 1  ■  [rj  1  [s  ]Lt] 

When  n>*2,  the  solution  obviously  becomes  indeterminate  since  the 
rank  of  the  matrix  is  then  less  than  the  order,  and  recourse  must 
be  had  to  least  squares. 


(15) 

(16) 


2.  The  Fundamental  Plane  Equations: 


Hie  general  approach  given  above  has  not  been  used 
for  computations  because  of  the  difficulty  of  solving  the  equations 
with  small  mechanical  computers.  Instead,  as  much  of  the  compu- 
tation  as  possible  is  done  on  the  fundamental  plane,  X1  x2  .  When 

the  equations  are  formulated  in  xi,  coordinates,  a  tremendous 
simplification  takes  place,  r*  disappears  formally,  since  r* 

is  parallel  to  x^;  projected  onto  xl  x2  becomes  identical  with 
the  observed  radius  of  the  moon,  and  rm  projected  onto  X1  x2  shrinks 
approximately  to  the  site  of  the  earth's  radius  or  less.  The  radius 
vector  from  the  spheroid  center  to  the  observer  is  also  projected 

onto  x  x  ,  and  the  combination  of  the  two  equations  (15)  and  (16) 
on  the  fundamental  plane  is  then  very  simple. 

Let  0*  be  the  projection  of  onto  the  plane  x*-  x^,  and 
let  £  ^  be  the  vector  along  the  x3-axis.  Then 


which  maps 


7 


-  <  z  ?>  v] 

(Vb  IH 
h 


(17) 


(18) 
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(19) 


The  simultaneous  solution  of  (14)  amd  (18)  la  given  by 

MTr>  [BJ  -  tqTQsi  -  K.i\. 


This  la  a  second-degree  aquatloa.  To  make  it  aaalar  to  work 
with,  It  can  be  linearised  by  taking  the  differential  of  both 
sides,  to  give  the  first-degree  equation  in  C&€][  and  3 

T, 


ex' A 


■  M  -  g.j  [MI  -  [^.f 


(20) 


From  here  on,  the  treatment  la  either  by  the  alngle-alte 
method  or  by  the  equal-limb-line  method.  The  former  will  be 
treated  flrat. 

C.  Single -Site  Method  of  Survey 

The  equation  (20)  above  la  so  simple  in  form  that  it 
gives  a  misleading  idea  of  the  number  of  important  quantities 
which  actually  influence  the  results.  In  the  equal-limb-line 
method  this  is  not  too  important,  since  moat  of  these  influencing 
factors  are  removed  by  the  control  observations ;  in  the  single¬ 
site  method  it  is  very  important  to  identify  every  influencing  factor 
and  to  evaluate  its  else  in  order  to  be  able  to  decide  whether  the 
factor  can  be  ignored  or  must  be  solved  for  along  with  the  geodetic 
unknowns.  To  emphasize  the  presence  of  these  factors,  equation  (20) 
is  rewritten  as : 

where 

[i‘]= 

The  [u^J  are  the  factors  whose  Influence  on  the  rtaults  is  to  be 
determined.  A  large  part  of  the  evaluation  is  of  a  quantitative 
nature  and  will  be  covered  in  Fart  II,  Error  Analysis.  At  this 
point,  only  a  general  discussion  of  the  particular  factors  involved 
is  needed. 
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The  most  obvious  factors  ere  the  geodetic  corrections: 

AE*  *  AE,  the  correction  to  the  eesting 
2 

AE  *  AN,  the  correction  to  the  northing  (22) 

3 

AE  s.  AH,  the  correction  to  the  height  ebove  the 
spheroid, 

which  enter  into  the  variables  ^a£\.  The  [^^ml  give,  of  course, 
the  fectors : 

Af  <23> 

in 

Ar 

m 

defining  the  corrections  to  the  position  of  the  center  of  the 
moon,  end 

A  b,  (24) 

m 

&P. 

which  define  corrections  involved  in  loceting  the  occulting 
f eeture . 


Implicit  in  the  equetion,  however,  and  not  literally 
present  ere  e  number  of  other  quantities.  First,  there  are 
the  quantities 


(25) 


which  were  used  to  define  the  space-fixed  directions  of  the 
x1 -system.  Next,  there  ere  the  quantities 

AX1  (26) 

mo 

which  enter  into  the  analysis  because  the  equation  was  set 
up  as  if  everything  were  referred  to  the  observer's  reference 
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I 

spheroid  and  datum  at  P,  whereas,  in  fact,  the  position 
XJof  the  moon  is  taken  from  the  Brown's  Tables  of  the 
Motion  of  the  Moon  (reference  4)(  or  from  an  ephemeris  based 
on  these  tables,  reference  6),  and  these  tables  refer  [x^ 
to  a  system  which  is  at  least  approximately  geocentric  — 
that  is,  which  has  its  origin  at  the  center  of  gravity  of 
the  earth.  The  vector  j  is  therefore  the  displace¬ 

ment  between  the  spheroldoclfttric  system  and  the  putative 
geocentric  system  origins. 

Another  quanity  which  must  be  considered  is  AT,  the 
error  in  the  recorded  time  of  occultation,  which  enters  not 
only  into  but  also,  through  ,  into  [£]  .  The 

time  error  is  of  two  kinds:  instrumental  and  astronomical. 

The  former  arises  from  errors  in  the  mechanical  equipment 
used  in  recording  and  reading  the  time  associated  with  the 
occultation,  while  the  latter  arises  principally  from  the 
difference  between  the  time  scale  in  which  the  observations 
are  carried  out  (U.T.)  and  the  "time"  scale  associated  with 
thefx,,,^  given  in  Brown's  Tables  (E.T.).  This  difference  in 
scales  arises,  as  is  well  known,  from  the  irregular  rotation 
of  the  earth  (which  rotation  measures  U.T.)  as  compared  with 
the  revolution  of  the  moon  about  the  earth  (which  measures 
E.T.),  and  from  the  failure  of  Brown's  tables  to  represent 
the  motion  of  the  moon.  The  former  effect  is  apparently  very 
much  greater  than  the  latter. 

In  addition,  there  is  a  quantity  £  4  f* m\  analogous  to 
the  which  represents  the  displacement  of  the  visible 

or  geocentric  center  of  the  moon,  which  is  what  is  observed, 
from  the  center  of  mass  of  the  moon,  which  is  given  as  [xmJ 
in  Brown's  tables.  Very  little  is  known  regarding  the  actual 
value  of  >  only  the  Ab,,,  component  (i.e.,  the  latitude 

component)  is  known  well  enough  to  be  given  a  value  in 
computations. 

The  distance  itself  contains  several  different  types  of 
error  which,  however,  cannot  be  solved  for  but  must  either  be 
eliminated  by  a  least-squares  adjustment,  or  allowed  fdr  in  the 
reduction  by  use  of  the  best  available  information.  The  most 
Important  of  these  is  the  quantity  hn,  which  is  the  height 
of  the  occulting  feature  above  the  lunar  mean  datum,  hm  Is 
different  for  each  occultation,  hence  cannot  be  found  Trom  the 
analysis;  it  can  be  estimated  to  within  200  meters  from  Hayn's 
charts  (reference  5)  or  from  the  charts  of  Dr.  Watts  of  the 
U.S.  Naval  Observatory  (reference  6).  Only  the  latter  charts 
can  be  considered  to  be  accurate  to  ±  200  steters,  although 
Hayn's  charts  are  contoured  at  0'.'2  or  300-meter  intervals. 
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When  the  list  of  factors  is  further  investigated,  it  is 
found  that  the  AS  +1  A  S  ,  Ar  can  be  broken 

down  still  farther.  The  £4  «*>  4$'*]  will,  lf*ke  h^,  vary  from 
occultation  to  occultation-.  In  addition,  because  of  the  nature 
of  occult'sf.ons,  they  cannot  be  separated  from  the  gross  values 
of  •  Hence  in  the  solution  they  will  be  lumped 

with  £Ac(m,  A«fm,  Ar]#Ths  corrections  £40,,  A  Arm]  to 
the  noon's  position  are  themselves  influenced  by  the  errors 
in  the  parameters  of  the  moon's  motion.  These  parameters  may 
be  variously  chosen;  the  paper  (reference  2)  by  Pamella  A. 
Henriks en  selects  the  following  correction  factors: 


Ar  5 

correction 

to 

the 

Ac (l»  ■ 

correction 

to 

the 

A  e  - 

correction 

to 

the 

Ai  - 

correction 

to 

the 

U  3 

correction 

to 

the 

A?  - 

correction 

to 

the 

longitude  of  perigee 
longitude  of  the  node 
eccentricity  of  the  moon's  orbit 
angle  between  lunar  and  stellar  axes 
obliquity  of  the  ecliptic 
position  of  the  equinox. 


(27) 


The  AT  could  also  be  included  in  the  group  of  factors  involved 
in  the  moon's  motion. 

Taking  all  the  above  correction  factors  as  a  set  Tau^  1 , 
the  equation  (21)  can  then  be  used  to  determine  the  fau^J 
if  sufficient  values  of  4<T  are  available.  If  there  are 
J  observations  (J  -  1  to  J,  J  >  K) ,  the  observation  equations  are 

-[y]rKH-k] 


(28) 


(29) 


In  order  for  the  least  squares  adjustment  to  give  reasonable 
values  for  C  au^J  ,  there  tnust  be  sufficient  observations  to 
provide  an  approximately  Gaussian  distribution  for  the  £r* 
and  the  h  ;  i.e,,  J>>K, 
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or,  In  Abbreviated  fora, 


l*v]  -[»] 


(31.1) 


Final 


W +  [AXJ  • 


The  [xj  are  related  to  the  [js]  and  the  [ax]  to  the  [ak]  by 
equations  3,  4,  5  and  the  partials  of  these  equations.  However, 
the 


d[?3  a[xJ 


are  zero,  so  that 


(31.2) 


Ihe  problem  cannot  be  considered  solved  by  this  equation  since 
it  is  the  position  [jcj  of  the  observer  with  respect  to  a  geo¬ 
centric  system  of  axes  which  has  been  found,  and  this  is  of  no 


value  by  itself. 

[*J 


What  is  wanted  is  the  relation  of 


B 


Final 


to  the 


of  some  standard  datum  and  spheroid  as,  for 
example.  North  American  Datum  and  the  International  spheroid,  or 
British  Datum  and  the  Airy  spheroid.  In  order  to  get  this 
relation,  a  similar  series  of  observations  would  have  to  be 
carried  out  independently  on  the  standard  reference  datum  and 
spheroid. 


Let  be  the  geocentric  coordinates  of  a  point  in 

a  reference  system  different  from  that  of  [Xy,  and  let  [x^j 
be  the  coordinates  of  this  same  point  referred  to  the  reference 
datum  and  spheroid  of  that  system.  Then 


M  '  W  +  ’ 

and  the  [AXS^|  are  found  as  before  from  the  equation 


(32) 

(33) 


and  equation  (31.2) 
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Since 


Au 


k 


n 


when  two  or  more  datum*  are  involved  it  i*  best  to  solve  all  the 
observation  equations  simultaneously,  except  for  the  geodetic 
variables , 


D.  Equal-Limb-Line  Method  of  Survey 


Return  now  to  the  equation  (20)  of  section  B, 

_._Trr 

(TACT 


•  [[A] 


M-  K2 


The  A  <T  can  be  separated  into  two  parts : 

A<r  S  AT  ac r 


(34) 


(35) 


^  d”g  is  that  part  of  the  error  in  A(p  which  arises  from  errors 

in  the  geodetic  coordinates  of  the  observer,  while  A  <f"  is  that 

m 

part  of  the  error  in  A0~  which  arises  from  errors  in  the  observed 

— >  > 

lunar  radius  (and  hence  includes  error  in  r  ,  P  ,  and  r.). 

If  the  occultation  is  observed  at  two  different  sites,  P^  and  P2» 

*  "I  - [U1 , -  (UJ '!>*],- 

^  1  -[1*1,  -  tfJ2]*[Ef<J2- 

First,  point  will  be  so  chosen  that 

^[t]2  -  0.  (37) 

This  assumption  will  be  more  fully  discussed  in  Part  II,  in  the 
section  on  geodetic  errors.  For  the  present,  it  will  be  justified 
by  the  rather  obvious  statement  that  any  datum  and  spheroid  can  be 
chosen  as  the  basic  datum  and  spheroid  to  which  all  positions  are  to 
be  referred,  in  which  case,  points  already  geodetically  connected  to 
that  datum  can  be  considered  as  having  ng  geodetic  error.  There  will 
actually  be  some  error  associated  with  r2g  and  hence  with  r^J^,  and 
we  could  write:  L 

t  -  ±T,  (37) 

2  02  2 

but  at  the  moment  it  is  enough  to  state  that  the  datum  and 
spheroid  on  which  P2  is  located  are  the  datum  and  spheroid  taken 
as  known.  Then  P2  can  usually  be  located  anywhere  within  a 
fairly  large  area  without  violating  the  assumption,  and: 
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‘’I A  •[[&, 

r»  *  {w* 


-  0s'1) 

.  cortK^  (39.2) 


NoCe  that  DP  -  and  t.43  2  411(1  k#nce  L^^jl  *nd  refer  to 

the  obsenrlng  sites *  coordinates, while  itris  the  corrections 
and  ££  cl-d  to  d*|u*  coordinates  which  are  sought  .  The 

relationship  sstween  [&4  1  and  1 3  .  is  not  easily  expressed  in 
the  -system,  where  the  rigidity  or  the  geodetic  net  requires  that 


[ae]" 

[*dd 

(40) 

Secondly,  P2  is 

chosen  in 

such  a  manner  that,  at  the  time  t^, 

^*2  “ 

(41) 

and  therefore 

^m2  * 

<r*. 

This  is  in' practice  an  extreaely  restrictive  condition.  Lack  of 
accurate  knowledge  of  the  shape  of  the  noon's  surface  makes  it 
imperative  that  the  occupations  as  observed  from  P~  and  P.  have 
(.JLU 2*  bm2^  48  c^°8e  48  possible  to  (  bml).  £m  and  b  are 
the  seleonographic  longitude  and  latitude  of  the  point  on  the  moon 
at  which  the  occultation  is  seen  to  occur.  When  the  moon's  llbratlons 
in  latitude  and  longitude  are  resolved  into  components  with  axes  paral¬ 
lel  to  and  x*,  it  is  seen  that  only  that  part  whose  axis  is  parallel 
to  x^  can  be  accounted  for.  This  component  is  a  rotation  in  position 
angle  Using  a  theorem  discovered  by  O'Keefe, 


dx 

dt 


2  TT 


sin  if 


(42) 


where  T„  is  the  moon's  sidereal  period;  if  is  the  Inclination  of 
the  moon's  axis  of  rotation  to  the  fundamental  plane;  and  "p  is 
the  poaltlon  angle  (reference  7).  Usd  of  this  equation  helps 
determine  the  coordinates  of  ?2,  since 


sin  ^ 


are  fixed, 


becomes  a  function 


(43) 


Once  the  coordinates 
of  the  tlste: 
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This  (second)  condition  therefore  limits  the  choice  of  points 
jwailable  to  ?2  to  a  line  which  is  the  locus  of  the  projection  of 
r^iC  fai,  t>m^)  onto  the  fundamental  plane. 

Once  ?2  has  been  chosen  to  satisfy  the  second  condition,  the 
final  equation  is  derived  by  subtracting  equation  (39.2)  from 
equation  (39.1).  Then 

^  tr  .1"  K]  i  ^  J  i 
!  -  2  “ 

L  <r. 
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i. 


(46) 


But  at  t^  and  tj,  beaause  of  the  choice  of  the  coordinates  of  Pj 


W 


Hence  as  can  also  be  seen  inraediately  from  the  diagram. 


»-» 

1 

H 

_ ft 

T 

1  >  1 

tfh) 1 

i_a _ , 

iT _ 1 

2  "  [Ami  2 

T_ 

a  i>  1 

L  *r 

and  therefore 

A  LN 

(47) 


»]r  M,1 

T 

At 

a^i  " 

<f 

1 

(48) 


This  is  the  basic  equation  for  the  equal-limb-line  method  of  survey 


giving  the  geodetic  corrections 


[t] 


as  a  function  of  the  difference 


between  the  observed  radii  of  the  moon  at  sites  and  P2. 

The  [a  must,  of  course,  be  related  to  the  which 

are  the  corrections  to  the  datum  of  P^  op  the  surface  of  the  earth. 
The  relations  are  provided  by  the  equations 

I 

.-l 


LC-  LhJ  LvHx]  +  [*  <J  Ll mJ  W  L 

L*»]  -  [*fj  LyJ  M+  Cvl 

Here  £xj  and  [^mj  refer  to  the  usual  astronomic  coordinate  system 
with  center  at' the  center  of  spheroid,  X^-axis  along  the  polar  axis, 


(49) 


are 


and  X*  parallel  to  the  meridiem  of  Greenwich.  j^E^j  and 
topocentric  systems  with  origin  at  the  observer's  local  reference 
point,  E^,  along  the  normal  at  that  point,  and  jjglj  tangent  to 

a  parallel  of  latitude.  The  'a  are  the  obvious  rotation 

matrices  which  rotate  from  the  £xj  and  systems  to  the 

fundamental  pldne  system  as  indicated  by  subscripts. 
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Writing 


M  ■  [*  5*]  [?  k) 

(50) 

M  -  [*  <fJ  [r/*J  W1  [r^j1 

and  differentiating, 

[A£]  -  0]  [ah]  ♦  [.]  [ax] 

(51) 

+[1¥Jr*]uU"J  +  [1^JW[Au“] 

where  £u*] 
that  order. 

■  X,  hj  T  for  n  -  1,  2,  3,  li,  5, 

in 

The  quantity  Au^  can  be  separated  into  two  parts: 

Au2  *  Au2*  +  Au^2 

6u  -  “f  - 
*  * 

(52) 

where  ©is  the  G.S.T.,  A  ©  is  certainly  less  than  0.01  or  5  meters, 
and  can  be  neglected  in  the  computations.  It  would,  in  any  case,  be 
random  in  character  and  would  therefore  be  among  the  quantities 
appearing  in  the  r.  m.  s.  error  after  the  least  squares  adjustment. 

The  errors  AX,  A(p ,  and  Ah  are  survey  corrections;  these  could 
presumably  be  found  if  enough  observations  were  available.  For 
most  sites,  however,  the  total  number  of  equal-llmb-lines  cannot 
be  expected  to  be  more  than  six  or  seven  in  number,  and  it  would  not 
be  desirable  to  include  AX  and  A (p  among  the  unknowns  to  be  found. 
Since  the  two  quantities  are  certainly  not  greater  than  A@,  these  also 
will  be  Ignored  in  the  solution.  This  leaves  the  unknowns  && *  and  A(X4 
and  here  again  an  analysis  of  the  expected  magnitudes  shows  that  the 
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error*  contributed  by  each  will  not  exceed  a  few  meters.  The 
total  error  caused  by  dropping  [au^]  altogether  --  i.e.,  assuming 


tu°]  “  [°] 

is  less  than  10  meters,  which  can  be  tolerated.  Equation  (51) 
therefore  becomes 

+  [a]  [axJ 


(53) 


(54) 


all  quantities  being  evaluated  at  point  Pi  (and  therefore  at  time  t^). 
The  quantities  u  and  both  represent  corrections  to  the 

position  of  P^;  the  latter,  however,  is  the  error  in  an  arbitrarily 
selected  quantity,  the  assigned  coordinates  of  the  origin  of  thp 
system  at  Pi,  and  hence 


(55.1) 


The  final  relationship  is  therefore 

[A  I,]  -  0J  [ag 

and  the  elements  of  are  the  direction  cosines 

cosCjS  E^) . 

When  this  is  substituted  back  into  the  basic  equation  (48),  there 
results : 


(55.2) 


A  IT, 


gl 


M.-fc), 


H  H-J 


or 


A<rgi  -  M*  03  i 


(56) 


(57) 


where 


is  a  vector. 
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If  J  observations  are  made  from  P  and  P-,  the  coordinate 
corrections  to  E^  are  given  by  the  equation: 


[a‘H 


The  variances 


associated  with  the 


are  given  by: 


(58) 


(59) 


where  the  S  is  the  variance  of  the  observation  error  distribution. 
These  [s1!  are  referred  to  the  E^-system  by  rotation  of  the  axes; 
the  off-9iagonal  variances  can  be  reduced  to  zero,  since  the  [s 1  j 
matrix  is  symmetrical.  Hence,  in  the  two-dimensional  case,  it  is 
possible  to  find  a  matrix  V-1  such  that 


,t  yi 

\  Vk* 


(60) 


The  rotation  ®  is  given  by 


where 


[$  ■  [ 

*]T[si]LRl 

r-i  B 

cos  ©  sin  @ 

LRJ 

-  sin  (h)  cos  (h ) 

(61.1) 


(61.2] 


or,  when  the  variance  ["AE^J  is  explicitly  included, 

cos  ®  sin  @  ol 

-  sin  ®  cos  ® 


(61.3 


0  0  1 


These  equations  define  the  semi -major  axis,  vj,  and  the  semi -minor 
axis,V2>  of  the  ellipse,  as  well  as  the  angle  lO  between  the  semi 
major  axis  and  the  E^-axis. 
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These  equatioms  are  easily  generalized  to  the  n-diaensional 
case,  which  is  covered  briefly  in'  Appeadix  E. 

Because  of  the  few  observations  available  from  amy  one 
site,  it  is  highly  desirable  to  restrict  the  solutioasat  first 
to  those  uakaowas  whose  errors  are  greatest.  The  coapoaeat  is 
heace  usually  set  arbitrarily  equal  to  zero: 

ae3  -  0, 

1  2 

aad  the  solutioa  Halted  to  the  determination  of  AE  aad  AE  , 
or,  what  is  the  saae  thing  for  the  purpose  of  computation,  A  A  cos 
and  A  .  in  this  connection  it  is  worth  pointing  out  again 
that  the  occultation  data  can  provide  no  information  about  the  x^~ 
component  of  events.  In  occultations  involving  point  P^  and  P2, 
which  are  2000  km  or  more  apart,  the  many  practical  restrictions 
on  the  observation  conditions  —  per  cent  illualnation  of  the  moon, 
time  after  evening  twilight  or  before  aiming  twilight,  etc.  -- 
often  combine  to  force  the  observations  at  P^  to  be  made  either 
near  the  horizon  or  near  the  zenith.  In  the  first  case  the  E^- 
component  will  be  poorly  determined  by  the  solution,  and  in  the 
second  case  the  E-* -component  would  be  poorly  determined. 

Up  to  now  it  has  been  assumed  that  the  corrections  are 

being  applied  directly  to  the  site  coordinates.  For  a  proper 

solution  the  corrections  should  be  applied  to  the  datum  ^  . 

coordinates,  since  it  is  assumed  in  the  derivation  of  the  equations 

that  the  site  coordinate  errors  result  from  the  datum  coordinate 

errors  (survey  errors  being  neglected;  as  explained  earlier). 

Where  the  observation  sites  are  only  a  hundred  kilometer*  apart, 

no  significant  error  is  Introduced  by  neglecting  the  difference 

between  6EiJobserver  «<*  &EiJdatum.  At  di8t"ces  much  greater 

than  this,  such  a  simplifying  assumption  cannot  be  aade,  and  a 

transformation  from  TeO  .  to  [e*Q.  __  i*  necessary. 

observer  *-  da  tun 


This  transformation  can  be  made  by  using  the  formulas  of 
Hrlstow  (reference  8)  connecting  the  change  in  position  of  a 
point  P  to  the  change  of  position  of  a  point  P2.  These  formulas 
can  be  put  into  the  form: 
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(61.4 


“  “ 

dA  cos  0 

d  A  cos  ^ 

df 

» 

df 

d  A 

d  A 

d  S 

- 

T 

d  S 

d  a 

d  a 

d  f 

-  - 

d  f 

—  - 

obs . 

- 

datum 


where  iTjis  a  function  of  the  separation  A0,  AA.  between  P-  and  P  . 
For  all  practical  purposes,  df  and  da,  the  changes  in  the  flattening  * 
and  in  the  semi -major  axis  of  the  spheroid,  can  be  set  equal  to 
zero  when  only  a  small  number  of  observations  are  available.  When 
many  observations  are  available,  covering  a  large  area  of  the  world, 
a  solution  including  df  and  da  could  be  significant.  Ihls  is  not 
the  case  at  present. 

Since  P-  is  tied  to  Pj  by  geodetic  survey,  it  may  be  assumed 
that  and  DAd-tum  are  not  to  be  changed  by  the  translation.  The 
resulting  equation  is  therefore: 


AE1  ~ 

ae1' 

AE2 

B 

X 

AE2 

(61.5) 

AP3 

observer 

L 

datum 

where  LT^Jis  the  matrix  |_TJ suitably  modified  (see  Appendix  F). 
This  transformation  can  now  be  Introduced  into  (57)  to  give: 


■  mVj  oi  d.tu.  i 


(61.6) 
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which  becomes 


ML.  •  {WT[*]  <«■»> 


analogously  to  (58). 

Th#  above  equations  57)  -  (61 . 7 )_/  have  been  forttulated 
for  the  case  where  only  two  points,  P  and  P  were  involved,  and 
where  Pj  was  the  "unknown"  point  and  P»  the  control  point.  In  an 
occultation  survey  prograa  covering  a  large  area,  the  distinction 
between  "unknown"  and  "control"  points  will  become  largely  artificial, 
with  a  particular  site  being  referred  sometimes  to  an  "unknown", 
sometimes  to  a  "control"  datum.  A  network  of  occultation  sites  will 
therefore  be  set  up  in  the  area  being  surveyed,  and  the  equation 
(61.7)  must  be  swdlfied  accordingly.  In  the  light  of  the  discussion 
in  Section  C  of  the  Single-Site  Method  Survey,  it  is  obvious  that 
an  adjustment  of  an  equal -limb -line  network  in  which  no  one  point 
is  held  fixed  (i.e.,  for  vhichfA^J*  0,  will  not  have  any  meaning.) 
Some  point  in  the  network  must  be  assumed  known,  so  that  for  it 


£k]  -  0  (62) 

The  occultation  equations  will  then  be  /referring  back  to  equation  (57}/: 

AS*  -  A*"  AE1 

J  ji  n 

i  (63) 

AE  a.  0 
k 


and  the  solution  will  be  given  by: 
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or  the  suitable  variation  according  to  (61.7).  The  variances 
are  then  found  as  before,  by  an  obvious  extension  of  the  previous 
equations . 

Suppose  D  distinct  datums  are  involved,  and  let  these  datums 

be  numbered  d  ■  1,  2,  3,  .  D.  Each  datum  is  to  be  fixed  on 

the  basis  of  observations  from  s  sites,  referred  to  that  datum, 
so  that 


Furthermore,  each  site  will  have  observed  q  occultatlons 

sd 


q  .  *1  ,2  ,3  , 

Msd  sd  sd  sd 

To  simplify  the  notation,  this  can  be  written: 


Then 


where 


d  «  1,  2,  3,  ....  D 

sd  *  Id,  2d,  3d . sd 

qsd  **  lsd,  2sd,  3sd . Qsd. 


(64.1) 


t,e  /  s,  d. 


From  (61.5)  and  (61.7)  and  the  above  equation  (64.1),  it  then 
follows  that  the  final  set  of  equations  can  be  written  as: 
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qsde 


qsd 


sd 


- 

S- 

m  "WE- 

a 

.  - 

+ 

B1 

qsd 

■r  41 

;  t1 

te 

H  X 

Ae1 

e 

»  -s 

(64.2) 


r  1 

"  * 

* 

i 

qsd 

i 

. 

I 

t4  *0 

W 

4 

+ 

lu 

i 

^E1 
_  e_ 

or,  enlarging  Che  vectors  and  matrices  and  Inverting, 

-{[  k]t[k]}'1  [Kf  [40i 

where  the  w- 


(64.3) 


vector  now  contains  all  the  datum  corrections. 


When  more  than  two  sites  lie  on  the  same  equal-limb-line, 
and  only  two  sites  are  involved,  it  is  obvious  that  the  total 
addition  to  the  geodetic  information  is  not  equal  to  the  number 
of  permutations  available,  but  is  actually  a  good  deal  less. 

These  cases  can  be  taken  care  of  by  weighting  the  observations. 

For  example,  the  most  common  case,  where  two  observations  are  made 
from  one  end  of  a  line  and  one  observation  from  the  other, 
could  lead  to  three  equations  involving. 


4<r 


12’ 


23’ 


a  <r 
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where  Pj  and  P-j  are  on  one  datum  and  Pj  on  another.  The  quantity 
6k<7"23  of  no  importance  to  the  adjustment  scheme;  ^  ^12  an<* 

A  <r“ 3i  are  of  importance  in  relating  the  two  datums  and  the  single 
observation  there  and  give  two  equations.  The  point  Pj  is,  however, 
common  to  both  equations.  Whereas,  if  the  equations  were  being 
solved  by  the  method  of  individual  occultations  (Part  II),  there 
would  be  three  equations  corresponding  to  the  three  values  A  O' 

&  0"g2t  ^°"g3*  48  opposed  to  the  four  equations  if  there  were8four 
independent  observations .  Hence  each  of  the  two  ^  ^ j  equations  is 
to  be  weighted  by  3/4  for  a  total  of  3/2.  A  table  for  relative 
weighting  of  observation  equations  referring  to  the  same  equal-limb- 
line  can  be  drawn  up  in  the  following  manner  to  take  care  of  the  most 
common  cases. 
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RELATIVE  WEIGHTING 


No.  of 

Sites  on  Line 

No.  of 
Equations 

Weight 

Total 

Weight 

Datum  1 

Datum  2 

per 

Equation 

1 

1 

1 

1 

1 

2 

2 

■ 

1.5 

1 

3 

3 

1  1 

1.5 

2 

2 

3 

0.75 

2.25 

written  as 


E .  Other  Forms  of  the  Occultation  Equations 

The  fundamental  relation  between  and  w  can  also  be 

>T|  sin'll 

(65) 

+  [Vf"jJT['os71'j  ^ 


m  -K  -Afe./L,iniy 


T 

(since  *2  =■  ~  -^U), 


Going  through  the  same  procedure  as  before,  it  is  easy  to  show  that 
the  final  equations  are,  formally, 


[aaj]  -  [asj]  ♦  OjfoTy 


(66) 


1  2 

where  £c»<3>Tj  has  the  usual  meaning.  fASjJis  the  x  x  -plane 
component  of  the  correction  to  ^E^,  and[Sj  ]T  is  t^e 

x*  x^-plane  component  perpendicular  thereto.  The  equations  can 

be  considered  as  giving  the  correction  components  [£E]or  [at] 
in  still  another  system,  the  s^-system,  where  the  axes  are  designated 
as 


(67) 


31 


The  e-axle  make*  an  angle  7"  with  the  *  -axis,  end  the 
x-axle  makes  an  angle  7  with  the  &  -axle.  Then 


[*0  ■  W  [*] 

«  [q]  [*] 


[■][■] 


(68) 


The  rotation  matrlcee £  rJ  ,  [rj  and  [qJ  then  create  vectora 


R,  r,  and  Q,  whoae  sum  la 


,  which  can  ltaalf  be  conaidered 


as  a  component  of  a  3-vector  M  whose  components  and  d~^ 

are  identically  zero.  The  rest  of  the  solution  can  be  carried  out 
just  as  before  and  is  no  different  mathematically  from  tha  previous 
methods.  The  interest  lies  in  the  fact  that  can  be  transformed 

into  a  correction  in  the  E^  E^-plane, 

T 


[ri]  M  ■  MH  +  MM  CM 


(69) 


dr 


AP 


where  £r  and  fcP  are  the  radial  and  tangential  components, 
respectively,  of  tha  corrections.  If  only  one  of  these  is  solved 
and  the  other  is  ignored,  then  the  question  is  whether  to  set 
Ar  or  dP  equal  tozzero  for  a  particular  occultation.  If 
dp  ■  0,  equation  (69)  defines  lines 

1  ■  (E  -  A  4r  sin  ^  )  sin  ©  (70) 

J  J  J  * 


and  the  corrected  position  of  the  point  P,  is  obtained  by  solving 
these  equations  for  the  intersection  of  the  straight  lines.  If, 
on  the  other  hand,  dr  ■  0,  equation  (69)  defines  circles  of  radius 
B,  AP,  and  the  corrected  E*  and  E^  lie  at  the  common  intersection 
of  these  circles. 

F.  Relative  Position  Method 


The  equations  and  methods  described  heretofore  are 
satisfactory  if  a  set  of  related  geodetic  positions  referred  to 
one  particular  datum  is  being  built  up.  This  is  especially  true 
when  the  set  is  global,  so  that  the  positions  esn  all  be  referred 
to  an  "absolute"  coordinate  system  (in  this  case  that  system  to 
which  the  lunar  orbit  is  referred).  Where  no  such  unified  set 
is  envisioned  or  possible,  however,  the  previous  equations  are 
open  to  the  objection  that  they  contain  Implicitly  the  assumption 
that  an  "absolute"  coordinate  system  is  involved  when  it  is  postulated 
that 
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for  one  of  a  pair  of  occultation  sites,  the  implication  is  present 
that  the  lunar  orbit  is  also  fully  defined  for  that  site.  This  is 
not  so,  and  the  resulting  error  will  not  be  removed  unless  a  global 
set  of  positions  is  involved. 

To  make  the  equations  useable  for  relative  position  deter¬ 
mination,  there  is  no  change  required  in  the  equations  other  than 
to  drop  the  condition 

AXj  -  0, 

and  to  substitute  for  it  some  relation  which  does  not  involve 
knowledge  of  the  absolute  coordinate  system.  The  required 
relation  is  arrived  at  through  the  obvious,  and  simply  proven, 
considerations  that  the  vectors  AX*  and  ^  Xg  must,  if  they 
represent  relative  shifts,  be  equal  in  magnitude  but  opposite 
in  signs j  the  unknown  relation  between  the  lunar  orbit  coordi¬ 
nate  system  and  the  local  geodetic  systems  is  the  same  at  both 
places.  The  resulting  equation  is  therefore: 

[aikl  +  [aik]  [Axi]  "  °21k> 


which  is  solved  for  the  [&11]  by  the  same  method  as  before. 

There  is  no  reason  why  this  equation  cannot  also  be  used 
in  setting  up  a  global,  connected  set  of  positions.  From  the 
relative  nature  of  the  Axi,  it  is  obvious,  of  course,  that 
the  end  product  will  be  "floating*  in  the  sense  that  two  (if 
is  involved)  or  three  parameters  must  be  arbitrarily 

specified  unless  the  connection  is  closed. 


The  preceding  discussion  has  covered  principally  the  data 
analysis  aspect  of  occultations.  For  the  prediction  of  occulta- 
tions,  the  same  equations  would,  of  course,  apply.  However, 
additional  work  must  be  done  to  select  the  stars,  compute  the 
"equal- limb-line  positions,  and  so  on.  The  procedures  have  been 
put  into  a  form  suitable  for  mechanical  calculation  by  the 
Department  of  Geodesy,  Army  Map  Service;  these  procedures 
will  be  published  in  a  forthcoming  report  entitled  "Mechanical 
Procedures  for  Occultation  Prediction." 
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PART  IL..  ERROR  ANALYSIS’- 


Introduction 

The  equations  which  govern  the  use  of  occultations 
for  surveying  are  incomplete  as  they  are  given  in  Part  I 
of  this  report,  since  they  do  not  state  the  range  of  errors 
which  can  affect  the  results.  These  errors  in  the  results 
arise  from  errors  in  the  constants  and  original  data,  and 
their  sizes  are  decided  by  both  the  sizes  of  the  input 
errors  and  the  equations  themselves.  An  analysis  of  the 
errors  therefore  concerns  itself  first  with  the  sources 
of  errors  in  the  input,  and  secondly  with  the  growth  of 
these  errors  as  they  pass  through  the  occultation  equation.;. 

As  is  seen  from  the  equations  in  Part  I,  the  sources 
of  error  are: 

1.  Star  positions 

2.  Moon  positions 

3.  Geodetic  coordinates  of  the  observer 

4.  Timing  of  the  occultation. 

Errors  in  star  position  and  moon  position  combine  with  the 
(unknown)  lunar  radius  at  the  point  of  occultation  .to  pro¬ 
duce  an  apparent  lunar  radius  which  is  greater  than  or  less 
than  the  standard  assumed  radius  of  the  moon.  Errors  in 
geodetic  position,  by  moving  the  observers  off  the  "equal - 
limb-llne"  positions,  affect  the  results  by  amounts  Ag  which  are 
approximately, 

Ag  - 

0  P 

where  JJJl  is  the  average  slope  in  the  vicinity  of  the 
dP 

occultation  feature  on  the  moon,  and  Aj£  is  the  error  in 
latitude.  Errors  in  timing  are  transformed  into  errors  in 
position  correction  which  are  approximately  equal  to  the 
time  error  multiplied  by  the  average  velocity  of  the  moon's 
projection  onto  the  observer's  horizon  plane. 
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II.  Error  Analysis 
A.  Error  Sources 

1.  Star  Position  Uncertainties: 

Hie  occultation  survey  method,  using  two  independent 
observers,  is  not  sensitive  to  reasonably  sized  errors  in  the  star 
coordinates.  Such  errors,  practically  inseparable  from  errors  in 
the  moon's  position,  cancel  out  almost  completely  in  the  first- 
order  approximation  to  the  exact  occultation  equations.  Sample 
computations  used  in  the  prediction  star,  positions  differing  by 
1",  have  given  paths  differing  by  only  1  or  2  meters  over  a  1000- 
kilcmeter  length.  However,  it  is  obviously  desirable  to  reduce 
all  errors  to  a  minimum  consistent  with  computational  ability  so 
as  to  produce  the  number  of  unknowns  that  must  be  accounted  for. 

The  source  of  the  star  positions  used  is  of  major  importance 
here.  Since  only  stars  occulted  by  the  moon  need  be  considered, 
attention  can  be  restricted  to  the  zodiacal  zone.  Furthermore,  it 
is  desirable  to  have  as  homogeneous  a  set  of  stars  as  possible. 

The  requirements  of  (a)  high  precision,  (b)  zodiacal  zone,  (c)  homo¬ 
geneity,  and  (d)  adequate  coverage  to  9th  magnitude,  effectively 
limit  the  choice  of  sets  to  (1)  the  Robertson  Zodiacal  Catalog 
(ref.  17)  and  (2)  the  Tale  Zone  Catalogs  between  +30°  and  -30°  (ref.  18) 
The  FKl*  contains  too  few  stars  to  be  useful  as  does  the  N30;  the 
Boss  General  Catalog  is  inaccurate  (one  second  of  arc  or  more  r.m.s. 
error).  Hence,  all  stars  used  have  been  from  either  the  Robertson 
or  Tale  Zone  catalogs,  with  corrections  being  made  where  necessary 
for  the  systematic  differences  between  the  two.  The  necessary 
corrections  are  made  by  using  the  conversion  tables  given  in  the 
Tale  Zone  and  Robertson  catalogs,  going  in  each  case  from  the  system 
in  question  to  the  Boss  General  Catalog  and  thence  to  the  appropriate 
second  system. 


Table  I.  Star  Position  Error 


Catalog 

Absolute  Value  of  p.  e. 

Tale  Zone  (1) 

<  055 

Zodiacal  (2) (3) 

<  055 

(1)  Evaluated  from  Tale  Zone  Catalogs  for  Zones  between 
♦30°  and  -30°  declination.  See  Table  II. 

(2)  Approximate  evaluation  from  50  random  samples  in 

catalog. 

(3)  Systematic  differences  between  the  Zodiacal  Catalog 
and  the  General  Catalog  run  from  about  +0Sl*5  to  -0922  in  right 
ascension  and  from  +0526  to  +0507  in  declination. 
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TABLE  II.  Prohibit  Error*  la  Diractioa 
Yalo  Zoao  Catalog* 


Zoaa 

oc 

S 

d«./dt 

d  SI  dt 

1960 

Total 

i 

-30  to  -27 

i  0.105 

±  0.105 

±  0.011 

t  0.010 

±  0.41 

-22,  to  -27 

±  0.105 

±  0.105 

±  0.011 

±  0.010 

±  0.41 

-20  to  -22 

±  0.105 

±  0.105 

-  ±  0;011 

±  0.009 

±  0.41 

-18  to  -20 

±  0.120 

±  0.105 

±  0.011 

±  0.009 

±  0.43 

-14  to  -18 

±  0.120 

±  0.105 

±  0.012 

±  0.009 

±  0.43 

-10  to  -14 

±  0.120 

±  0.105 

±  0.011 

±  0.011 

±  0.46 

-  6  to  -10 

1  0.105 

±  0.105 

±  0.008 

±  0.008 

±  0.32 

-  2  to  -  6 

±  0.105 

±  0.105 

±  0.008 

±  0.32 

+  1  to  t  2 

±  0.115 

±  0.115 

±  0.006 

±  0.007 

±  0.24 

+  5  to  +  \ 

±  0.115 

±  0.115 

±  0.007 

±  0.007 

±  0.27 

+10  to  +  5 

±  0.115 

±  0.115 

±  0.008 

±  0.007 

±  0.27 

+15  to  +10 

±  0.115 

±  0.115 

±  0.008 

±  0.006 

±  0.26 

+20  to  +15 

±  0.115 

±  0.115 

±  0.005 

±  0.006 

±  0.23 

+25  to  +20 

±  0.095 

±0.095 

±  0.005 

±0.005 

±  0.26 

+30  to  +25 

1  t 

±  0.095 

±  0.095 

±  0.007 

±  0.005 

±0.31 
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From  an  analysis  of  the  errors  introduced  by  each  of 
the  two  catalogs,  as  summarized  in  Tables  I  and  n,  it 
appears  that  there  is  no  significant  difference  between  the 
errors  of  the  two  catalogs  (proper  motions  being  taken  into 
account).  Except  for  systematic  differences  between  the  two 
catalogs,  there  is  no  hindrance  to  the  simultaneous  use  of 
both  catalogs  in  the  reduction  of  occultations.  The  different 
backgrounds  of  the  two  catalogs  (the  Robertson  Catalog  is 
"fundamental"  and  relates  to  an  earth-fixed  system,  while  the 
Tale  Zone  Catalog  is  differential  and  hence  more  "absolute" 
in  nature)  have  no  significance  as  far  as  the  use  of  the 
catalogs  for  occultations  is  concerned. 

In  summary,  then,  it  can  be  assumed  that  the  probable 
errors  in  right  ascension  and  declination  of  the  directions 
of  Ihe  stare  used  are  less  than  0?5>  and  hence  that  the 
standard  deviations  are  less  than  057$.  The  corresponding 
error  introduced  into  the  computed  geodetic  position  is  less 
than  2  meters. 

2.  Moon  Position  Uncertainties 


The  phrase  "center  of  the  moon"  is  inadequate  for 
defining  a  point  since  the  definition  used  is  operational  and 
actually  varies  with  the  method  being  studied.  There  are 
three  kinds  of  lunar  centers  which  are  of  interest: 

a.  The  center  of  mass,  or,  more  precisely,  the 
center  of  gravity. 

b.  The  optical  center. 

c.  The  operational  center. 

The  center  of  gravity  of  the  moon  is  at  present  defined  by 
Brown’s  equations  for  the  motion  of  the  moon  (ref.  U)  as 
modified  in  the  Improved  Lunar  Ephemeris  (ref.  19).  It  is 
a  set  of  coordinates  (latitude,  longitude,  and  parallax)  to  which 
the  center  of  the  moon  is  said  to  belong,  and  which  are  computable 
as  a  function  of  time  alone.  These  coordinates,  established  by 
the  equations  of  motion,  do  not  correspond  to  any  point  that  can 
be  located  by  observation,  because  (a)  the  theory  is  based  on 
values  of  basic  parameters,  such  as  the  size  and  shape  of  the 
earth,  whose  errors  are  not  known;  and  (b)  the  theory  refers 
its  predictions  to  the  gravitational  center  of  the  earth 
approximately,  while  actual  observations  do  not  have  access 
to  this  center. 
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The  right  ascension  and  declination  of  tha  noon  are 
given  in  the  Improved  Lunar  Ephcmeris  (ref.  19),  and  in 
national  ephemerides  of  1960  and  following,  to: 

07015  right  ascension 

O'.'Ol  declination. 

These  were  obtained,  however,  by  interpolation  in  the  computed 
longitude  and  latitude  of  the  moon  as  given  by  Brown's  equations, 
and  are  really  lisd.ted  in  precision  only  by  the  abilities  of  the 
computing  machine  used  to  evaluate  the  equations,  and  could  be 
carried  out  to  any  number  of  figures.  The  constants  entering 
into  the  computation  are  good  to  10~5  seconds  of  arc,  and  this 
together  with  neglect  of  terms,  machine  rounding-off  error,  etc., 
limits  computed  positional  precision  to  about  ±  07001.  As  is 
discussed  further  under  Time  Uncertainties,  the  precision  is 
here  also  the  accuracy.  The  reason  for  this  is  that  at  the 
present  time  all  discrepancies  between  the  computed  coordinates 
and  coordinates  obtained  from  observations  are  assumed  to  be  . 
caused  by  irregularities  in  the  observations,  and  primarily  in 
the  real  time.  As  far  as  occultatlons  and  most  other  obser¬ 
vational  techniques  are  concerned,  the  positional  uncertainty 
in  the  gravitational  center  of  the  moon  can  be  taken  as : 


(*> 

From  published  ephemerides 

<  0702 

(b) 

From  special  ephemerides 

<  07002 

(c) 

From  specific  computations 

07 

b.  The  only  center  for  the  moon  which  is  "directly" 
observable  is  the  optical  center.  Even  the  gravitational  center 
of  the  moon  is  not  free  of  optical  influences,  since  the  constants 
entering  into  the  equations  of  motion  were  determined  from  optical 
observations.,  It  is  true  that  the  Connection  is  tenuous*  since 
a  very  large  number  of  optical  observations  enter  into  the  con¬ 
stant  determinations,  so  that  random  variations  could  be  expected 
to  be  averaged  out.  Systematic  effects  do  appear  and  must  be 
accounted  for;  the  initial  work  by  Brown  (ref.  1,  4)  gives  a 
value  of  -075  for  the  difference  in  latitude  between  the  optical 
center  and  the  gravitational  center.  This  value  accounts  explicitly 
for  all  residuals  in  the  latitude  equation  not  otherwise  accounted 
for,  and  hence  is  not  entirely  satisfactory.  Furthermore,  there  is 
an  asymmetry  between  the  eastern  and  western  portions  of  the  lunar 
limb,  as  is  shown  by  existing  profile  measurements.  This  asymnetry 
does  not  show  up  in  the  equations  of  motion. 
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A  more  direct  way  of  finding  the  optical  center  is  to  fit 
circular  or  elliptical  arcs  to  lunar  profiles  and  to  resect  from 
these.  Or,  what  is  the  same  thing,  to  fit  a  circle  or  ellipse  to 
profile  measurements  to  give  best  fit  in  the  least  squares  sense. 
The  latest  and  probably  best  measurement  of  the  relation  between 
center  of  figure  and  center  of  gravity  is  that  of  Watts  and  Scott 
(ref.  23),  which  gives  -  0560  +  0508  for  the  difference  in 
latitude. 

The  value  given  in  the  American  Ephen^ris  and  Nautical 
Almanac  to  be  used  for  the  calculation  of  eclipses  is  -  055* 

c.  The  definition  of  the  moon's  center  is  that  in  which 
the  observations  (occultations,  eclipses,  etc.)  themselves 
provide  the  center.  In  this  way,  centers  peculiar  to  each  method 
arise.  For  example,  a  careful  statistical  analysis  in  193^ 
by  Mr.  G.  Reuning  of  some  2000  visually  observed  occultations 
gave  an  optical  center  minus  ephemeris  center  value  of  OS  ■  05U 
(ref.  2L).  Similarly,  eclipse  reductions  provide  temporary 
centers  pertaining' to  the  eclipse. 


TABIE  III.  lunar  Position  Uncertainties 


Uncertainty  in: 

Absolute 

Relative 

Difference  between  Optical  and 
Gravitational  Centers 

051 

~  0" 

Orbit 

(2) 

<  0502 

—  0" 

Right  Ascension 

<  05013 

~  0" 

Declination 

<  0501 

~  0" 

Radius 

(3) 

0 

0 

Height 

(U) 

0512 

«  0509 

To  compare  a  theoretical  position  of  the  moon,  as  defined 
by  the  Improved  Lunar  Ephemeris,  with  the  actual  position  of  the 
moon  as  defined  by,  say,  transit  or  moon  camera  observations,  is 
not  very  easy,  since  instrumental  errors  must  be  also  accounted 
for,  A  simpler  and  almost  as  valid  method  of  comparison  is  to 
compute  the  difference  between  the  Improved  Lunar  Ephemeris 
position  and  a  position  based  on  the  same  equations  but  with  more 
"modern"  values  for  some  of  the  fundamehtal  constants. 
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If  we  take  the  following  values  as  possible  values, 


a  6,378,150  meters 

f  1/298.3 

g£  978.036  cm/sec2 

M/fc  1/81.32 


for  the  semi-major  axis  of  the  earth-fitting  spheroid,  the 
flattening,  the  mean  gravitational  acceleration  at  the 
equator,  and  the  ratio  of  lunar  to  terrene  mass,  then  the 
moon's  horizontal  equatorial  parallax  as  given  by  Brawn's 
equations  requires  correction  by  a  factor 

1  -  6.72  x  10“5. 

while  the  lunar  ephemeris  computed  by  the  various  observations 
from  these  equations  requires  a  correcting  factor: 

1  -  2.05  x  10“5. 

Furthermore,  all  terms  with  the  coefficient  which  contain 
101  and  a,  must  be  multiplied  by  the  correction  factor 

1  -  1.805  x  10-3. 

A  list  of  those  terms  which  are  sufficiently  changed  to  affect 
the  position  by  05001  is  given  in  Table  IV.  Further  changes, 
in  addition  to  those  in  which  the  dynamical  parallax  occurs 
directly,  are  required  by  a  change  from  Brown's  value  of  1/29U 
for  the  flattening  to  1/298.3.  These  changes  are  of  the  size 

1  -  2.92517  x  10"2. 

Terms  affected  to  the  extent  of  05001  are  listed  in  Table  V. 
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TABLE  IV .  Llat  of  Terms  Affected  by  Changes  in  Major 
and  Minor  Axes  Dimensions  of  Ellipsoid 


Number 

Value 

Change 

Affected 

21 

- 

125.154 

+  0.250 

Longitude 

55 

- 

8.466 

+  0.017 

Longitude 

56 

+ 

18.609 

-  0.037 

Longitude 

57 

+ 

3.215 

-  0.006 

Longitude 

61 

+ 

18.023 

-  0.036 

Longitude 

62 

+ 

0.560 

-  0.001 

Longitude 

118 

- 

0.986 

+  0.001 

Longi tude 

119 

+ 

1.750 

-  0.003 

Longl tude 

120 

+ 

1.225 

-  0.002 

Longitude 

123 

+ 

1.267 

-  0.002 

i 

Longitude 

1 

129 

- 

1.089 

+  0.002 

Longitude 

138 

+ 

0.584 

-  0.001 

Longitude 

627 

“ 

0.978 

+  0.002 

Parallax 
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TAB  IE  V.  List  of  Terms  Affected  by  Change  of  Ratio 
of  Major  to  Minor  Axis  or  tillipaoid 


Number 

Change 

Term 

Affected 

1375 

1.  Arguments 

-  0.212 

L 

1376 

-  0.008 

L 

11*01 

♦  0.061 

'50' 

11*07 

-  2.807 

«/l 

Q fl 

11*08 

-  0.1*56 

U*09 

-  0.05U 

Ji 

2.  Periodic  Terms 

3 

-  0.007 

Longitude 

7 

-  0.025 

Longitude 

8 

♦  0.019 

Longitude 

25 

♦  0.002 

Longitude 

26 

♦  0.001 

Longitude 

33 

♦  0.001 

longitude 

51 

-  0.009 

Longitude 

31 1* 

-  0.007 

Latitude 

326 

♦  0.025 

Latitude 

328 

+  0.019 

Latitude 

339 

-  0.002 

Latitude 

.  The  Radius  of  the  Moon 


Aa  stated  previously,  the  principal  sources  of 
information  on  lunar  heights  are  the  charts  of  Hayn,  Weimer, 
and  Watts.  Using  these  charts,  it  becomes  a  matter  of 
indifference  what  the  mean  surface  is,  as  long  as  the  proper 
value  is  used  with  each  chart.  From  the  method  of  construction 
of  Watts'  profile,  the  major  error  is  introduced  in  the  process 
of  finding  the  mean  surface.  The  various  errors  can  be 
summarized  as : 

actual  surface  to  mean  surface  ±  0V07 

mean  surface  to  center  (optical)  ±  0.2 

optical  center  to  mass  center  t  0.1 

Total  ±  0723 

Although  the  computations  for  equal-limb-line 
occultations  presume  that  the  occulting  feature  is  located 
precisely,  the  actual  situation  is  quite  different.  A  major 
uncertainty  arises  through  the  lack  of  ability  to  account  for 
(a)  the  uncertainty  in  profile  caused  by  uncertainty  in  ground 
position,  and  (b)  the  uncertainty  in  profile  on  the  moon's 
edge  because  of  libration  in  longitude  between  occultations  at 
the  two  or  more  ground  stations. 

The  difference  between  the  assumed  position  of  a  ground 
point  and  its  position  referred  to  the  lunar  orbit  reference  system 
will,  in  all  cases  except  where  initial  survey  was  grossly  inadequate, 
be  less  than  3  km. ,  measured  along  the  reference  ellipsoid.  The 
greatest  uncertainty  will  occur  when  this  distance  is  along  the  lunar 
edge  and  when  the  slope  of  the  profile  in  that  region  is  maximal. 

The  radial  uncertainty  would  then  be 

A  r  *  ±  1500  meters 

max 

with  respect  to  the  radius  at  the  other  ground  sites.  Since 
statian  locations  are  changed  from  occultation  to  occultation  to 
diminish  the  difference  between  assumed  "true"  positions,  after 
one  or  two  occultations  the  radial  uncertainty  would  be  reduced 
to  correspond  approximately  to  the  uncertainty  in  the  occultation 
methdd  itself.  Furthermore,  occultations  are  selected  which 
occur  at  maria  or  similarly  flat  areas  in  which  the  maximum 
slope  is  closer  to  ±  5*.  A  realistic  estimate  of  the  radial 
uncertainty  maximum'  arising  from  station  position  uncertainty 
would  then  be 

Ar  >30  meters, 
max 
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Since  the  maximum  time  between  occultatioas  on  an  equal-limb  survey  is 
less  than  2  hours,  the  radial  uncertainty  le  the  sua  of  the  uncertainties 
from  the  longitudinal  libratioa  position  at  the  different  occultatlon  times 
In  2  hours  the  longitudinal  libratioa  is  less  than  9',  and  the  aaxlaua 
radial  uncertainty  attributable  to  introduction  of  a  new  feature  by 
longitudinal  llbration  would  be  between  i  15  meters  certainty,  and  aost 
of  the  tiae,  therefore,  much  less. 

5.  Observer's  Position  Uncertainties 


In  theory,  the  latitude,  longitude,  and  height  of  an  observer 
are  clearly  defined  quantities.  In  the  simplest  case,  a  spheroid 
(ellipsoid  with  two  equal  axes)  of  size  and  shape  approximately  that  of 
the  earth  is  selected.  A  point  on  this  spheroid  is  chosen  to  represent 
a  given  latitude  and  longitude  on  the  earth,  the  position  usually  being 
chosen  so  that  its  geodetic  coordinates  coincide  with  or  are  close  to 
the  actually  measured  coordinates  of  the  physical  point  on  the  earth. 

The  spheroid  is  then  oriented  (in  theory)  until  its  minor  axis  ispAtallel 
to  the  earth's  mean  axis  of  rotation  and  its  major  axis,  labelled  X', 
points  parallel  to  the  plane  of  the  Greenwich  meridian.  Again, in  theory, 
the  spheroid  is  translated  bodily  until  the  chosen  point  on  it  lies  on 
the  vertical  through  the  selected  point  on  the  earth,  and  it  is  then 
further  translated,  always  without  rotation,  until  the  separation  of  the 
two  points  is  equal  to  some  pre-selected  distance.  The  point  on  the 
earth's  surface  is  referred  to  as  the  datum  point,  and  its  latitude, 
longitude,  and  height  above  or  below  the  spheroid,  together  with  the  size 
and  shape  of  the  spheroid,  are  sufficient  to  define  the  geodetic  system. 
Subsequent  points  located  on  the  earth's  surface,  with  respect  to  the 
datum  point,  are  transferred  to  the  corresponding  points  on  the  spheroid. 

In  the  case  of  a  large  area  such  as  North  America,  the  process  of 
transferring  points  from  the  earth's  surface  to  the  spheroid  surface 
can  introduce  large  errors.  These  errors  arise  because,  at  great 
distances  from  the  datum  point,  the  tri angulation  errors  which  Increase 
with  distance  become  larger.  Furthermore,  the  distance  between  the 
surface  of  the  earth,  on  which  the  measurements  are  made,  and  the 
surface  of  the  spheroid,  to  which  positions  are  referred,  becomes  less 
well  known  with  increasing  distance  from  the  datum  point,  at  which 
the  spheroid-earth  difference  is  defined. 

Where  there  is  no  continuous  chain  of  survey  from  the  datum  point 
to  a  given  point  on  the  earth's  surface,  another  source  of  error  enters. 

This  is  the  unknown  amount  of  translation  needed  to  bring  the  spherdid 
proper  to  the  first  system  into  coincidence  (at  the  origin)  with  that  proper 
to  the  second  point. 
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The  majority  of  occultations  observed  so  far  by  the  equal- 
limb-line  method  have  been  such  that  one  of  the  two  observation 
points  lies  within  10  km  of  the  datum  point.  However,  a  sub¬ 
stantial  number  have  the  observation  points  at  considerable 
distances  from  the  datum  point,  and  hence  this  distance  must 
be  considered. 

Where  extensive  triangulation  networks  exist,  the  survey 
is  of  high-order  accuracy,  and  an  error  of  1  cra/km  in  the  distance 
between  points  on  the  network  itself  can  be  assumed.  Points  not 
on  the  network  are  located  to  within  ♦  1  meter  horizontally,  with 
respect  to  the  nearest  points  in  the  network.  For  a  point  1000  km 
from  the  datum  point,  the  error  will  therefore  be  approximately 
♦  10  meters,  horizontally,  the  local  error  being  insignificant. 

This  error  is  usually  that  in  distance  along  a  sea-level  surface. 
Because  of  the  not-very-well-known  separation  between  the  sea- 
level  surface  and  the  spheroid,  an  additional  amount  must  be  added 
to  show  this.  There  are  no  universally  valid  data  on  the  error  thus 
introduced;  there  are,  however,  excellent  approximations  available, 
based  on  assumptions  regarding  the  size  and  shape  of  the  sea-level 
surface,  which  make  a  value  of  +  meters  in  1000  meters  a  reasonable 
upper  bound. 

Although  even  a  very  long  triangulation  chain  represents  a 
high  order  of  accuracy  as  far  as  distance  is  concerned,  the  accuracy 
of  the  coordinates  at  the  end  of  the  chain  would  be  very  low  unless 
independent  means  were  available  for  controlling  the  direction  of  the 
chain.  This  control,  generally  afforded  by  other  chains  orthogonal 
to  the  first,  and  by  periodic  determinations  of  astronomic  azimuth 
along  the  chain,  can  be  assumed  for  the  purposes  of  the  analysis 
to  be  present.  It  may  be  noted  that,  in  almost  all  cases  where 
occultations  have  been  observed,  the  observation  sites  have  been 
well  within  1000  km  of  the  datum  point  and  usually  within  100  km 
of  the  datum  point. 

These  errors  or  uncertainties,  together  with  a  number  of 
others  of  ouch  less  importance  arising  from  errors  in  the  geoid 
above  spheroid  and  other  values,  are  all  considered  as  lying 
within  one  datum  system.  Another  large  error  arises  from  the 
unknown  relationship  between  positions  measured  on  spheroids 
which  are  located  in  different,  unrelated  datum  systems.  Because 
of  the  definition  of  the  spheroid  and  datum  or  because  of  the 
mathematical  procedures  used  in  computation  of  coordinates,  all 
spheroids  are  similarly  oriented.  The  errors  therefore  arise 
from  the  errors  in  measurement  of  the  distance  between  spheroid 
centers.  There  is,  of  course,  no  value  known  for  the  error  in 
this  distance  in  any  case  of  interest;  only  the  precision  is  known. 


Furthermore,  since  occultation  surveys  depend  on  e  knowledge  of 
the  luner  orbit  to  some  extent,  the  error  in  the  distance  between 
the  spheroid  center  end  gravity  center  is  of  importance.  A  value 
of  ±  200  meters  is  adopted  here  for  the  quantity.  This  value  is 
probably  greater  than  necessary,  but  not  unreasonably  so.  By  the 
very  mature  of  occultation  surveying,  this  particular  error 
appears  entirely  in  the  unknown  site  coordinates,  which  have  to 
be  corrected  in  any  case.  The  effect  of  the  error,  together  with 
the  errors  within  the  datum,  is  to  cause  a  displacement  of  the 
unknown  site  along  the  lunar  limb  and  perpendicular  thereto, 
which  Mans  that  the  unknown  site  is  going  to  obtain  a  final 
correction  whose  error  is  proportional  to  the  amount  by  which  the 
actual  feature  observed  on  the  moon  differed  from  that  predicted 
by  the  equal-llmb-line  calculations.  This  was  discussed  earlier 
in  the  section  on  lunar- position  uncertainties. 

6*  Time  Uncertainties. 


The  principal  sources  of  tlM  uncertainties  are  listed 
in  Table  VI,  together  with  estimates  of  the  uncertainty  magnitudes. 
The  uncertainties  arise  from: 

a.  The  nature  of  time  Itself  and  its  definition. 

b.  The  determination  of  tiM  and  its  propagation 
to  the  receiver. 

c.  The  reception,  amplification,  and  recording  of  tine. 

d.  The  propagation  of  the  light  signal. 

e.  Reception,  amplification,  and  recording. 

f.  Measurement  of  the  record. 

Total  time  uncertainty  is  found  to  be  less  than  0?01.  In  spite 
of  the  many  still  unanswered  questions  in  a  check  of  the  above 
list,  this  value  of  OfOl  is  considered  reliable  and  probably  close 
to  an  upper  bound.  The  reason  is  that  the  largest  part  of  this 
uncertainty  arises  from  the  Masuresmnt  resolution  limit,  which 
was  taken  as  0.5  am  upper  bound,  and  which  is  accessible  to  veri¬ 
fication,  and  from  the  radio  propagation  time  uncertainty  of  0f002, 
which  also  can  be  evaluated,  although  with  slightly  less  confidence. 
The  unknowns  in  the  tiM  problem  are  at  least  an  order  of  magnitude 
SMller  than  this. 
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TABLE  VI .  Time  Uncertainty  Estimates 


Absolute 
Delay  _(  J US) 

Relative 
Uncertainty  (iL. 

1. 

Delay  in  Electronic  Circuitry 

(1) 

a. 

Photomultiplier 

(2) 

0.1 

0.004 

b. 

Cable,  photomultiplier 

7000 

2 

to  amplifier 

c. 

Amplifier 

d. 

Cable,  amplifier-recorder 

<500 

<  10 

e. 

Recorder 

f. 

Radio  recevier 

<  500 

<.  10 

g- 

Filter,  1000-cps  band  pass 

2. 

Delay  in  Mechanical  Components 

(1) 

a. 

Recorder  pen  lengths 

1000 

b. 

Galvanometer  inertia 

(3) 

10,000 

<  1000 

c. 

Chart  speed 

<  1000 

<  1000 

3. 

Optical  Uncertainties 

»• 

Light  travel  time 

1,000,000 

■*  0.1 

b. 

Diffraction  edge  position 

(4) 

700 

<  100 

4. 

Radio  Propagation  Uncertainties 

Radio  propagation  time 

30,000 

<  2000 

5. 

Direct  Time  Errors 

a. 

Emission  time  correction 

(5) 

50 

<  1 

b. 

Ephemeris  Time  correction 

(6) 

10,000 

1 

c. 

Polar  motion:  rotation 

of  earth 

(5) 

500 

<■  5600 

6. 

Measurement  Errors 

a. 

Resolution  (0.5  mm) 

4000 

1 

b. 

Cross-chart  correction 

(0 . 1  mm) 

800 

1100 

TOTAL 

<8400 
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(1)  Tests  have  been  made  to  measure  the  time  delays  present 
in  the  equipment  between  signal  (light  or  radio)  reception  and  recording. 
The  system  paths  considered  are  shown  in  the  following  diagram: 


f SIGNAL  Pi 


_ ^ 

k _ 

RADIO 

1000  CPS 

FILTER 

PHOTOCELL 
- 7f - 


AMPLIFIER 


AMPLIFIER 


SIGNAL  IN 


■^RECORDER 
PART  I 


RECORDER 
PART  II 


RECORDER 
UNIT 


The  amplifiers  are  of  the  same  type;  the  recorder  units  (galvanometer  and 
pens)  are  of  the  same  type,  housed  together,  and  with  common  chart  drive 
and  chart.  Figures  5  and  6  show  the  two  test  methods.  Variants  of  these 
methods  were  also  used,  with  individual  components  of  the  systems  being 
tested  as  well. 

(2)  The  transit  time  in  nine-  or  ten-stage  photomultiplier  tube 

is  around  10-?  to  10-  seconds,  but  the  variation  in  this  is  many  orders  of 
magnitude  less.  Engstrom  (ref.  28)  and  Karoiuf  (rer  ? ) ,  among  others, 
give  values  for  the  transit  times.  There  is  some  caution  necessary  in 
the  evaluation  of  transit  times.  A  single  electron  at  the  input  to  an 
eleven?stage  tube  will,  according  to  the  literature,  produce  a  pulse 
of  10-  -second  width  at  half  amplitude.  The  signal  transit  time  is  the 
time  cf  start  of  the  pulse,  but  the  measured  time  will  depend  upon  the 
sensitivity,  averaging  properties,  etc.',  of"" the  measuring  instrument,  and 
could  be  the  pulse  half -width  point. 

(3)  The  galvanometer  response  time  was  evaluated  separately 
and  as  part  of  the  over-all  system.  Appendix  A  gives  the  derivation 
of  the  pen-reaponae  equations. 

(b)  The  moon  acts  as  a  diffracting  plane  to  produce  a 
diffraction  pattern.  The  effect  of  this  on  the  shape  of  the  occulta- 
tion  curve  and  on  the  time  determination  is  discussed  in  Appendix  B. 

(5)  The  radio  signals  are,  for  obvious  technical  reasons, 
seldomly  emitted  at  the  instant  which  they  are  supposed  to  indicate. 

The  actual  times  of  emission  are  determined  by  monitoring  of  station 
broadcasts  by  time  centers  such  as  Greenwich  Observatory;  the 
necessary  corrections  are  published  periodically  in  time  bulletins. 

The  major  source  for  such  corrections  has  been  the  U.  S.  Naval 
Observatory  Time  Bulletins.  The  bulletin  values  must,  of  course. 
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TEST  SET-UP  FOR  MEASURING 
DELAY  IN  TIME-SIGNAL  CIRCUIT 


SYSTEM  BEING  TESTED 


FIGURE  5 


Light 

source 


-case 


TEST  SET-UP  FOR  MEASURING 
DELAY  IN  STARLIGHT  CIRCUIT 


PHOTOMULTIPLIER 


FIGURE  6 
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be  corrected  for  time  of  transmission  from  the  station  to  the 
observatory}  since  only  the  difference  in  time  of  reception  by 
two  equal- limb -line  time  observers  is  involved,  this  transit¬ 
time  correction  is  of  no  importance  where  the  same  radio-station 
signals  are  used. 

Further  corrections  must  be  made  for  the  variation  of 
latitude  caused  by  the  motion  of  the  pole  and  for  the  seasonal 
variation  in  the  rotation  of  the  earth.  The  necessary  correction 
constants  are  computed  by  the  Bureau  International  de  l'Heure 
(ref.  27)}  preliminary  values  are  provided  in  the  time  bulletins. 

(6)  The  observations  of  occultations  are  related  to  tine 
through  radio  signals  sent  out  by  radio  stations ,  or  to  locally 
generated  time  signals  from  a  clock  which  is  calibrated  against 
the  radio  signals.  The  problem  of  accounting  for  the  difference 
between  the  time  recorded  by  the  occultation  equipment  and  the 
time  -  Universal  Time  -  at  which  the  time  signal  should  have  been 
emitted  is  complex,  but  the  individual  steps  necessary  to  the 
solution  are  well  understood.  This  is  not  the  case  with  the  relation 
between  Ephemeris  Time  and  Universal  Time.  Universal  Time  and  its 
variants  Universal  Time-1  and  Universal  Time-2  lend  themselves  to 
precise  operational  definition;  Ephemeris  Time,  on  the  other  hand, 
has  no  physical  existence,  is  defined  theoretically  in  terms  of  an 
unrecoverable  unit,  and  can  be  related  to  other  time  standards 
only  through  empirical  methods.  There  is  no  theoretical  relation¬ 
ship,  for  instance,  between  the  Ephemeris  Time  second  and  the 
second  defined  by  an  atomic  molecular  oscillator. 

The  unit  of  Ephemeris  Time  (ref.  19)  is  the  tropical  year 
defined  by  the  mean  motion  of  the  sun  in  longitude  at  January  0, 

1900,  12“  E.T.  The  ephemeris  second  is  defined  as  (31,  556,925.97U7) 
times  the  topical  year  at  the  above  epoch.  Since  Ephemeris  Time 
is  designed  explicitly  to  give  positions  of  the  major  solar  system 
bodies  which  agree  with  observations,  the  relation  of  Ephemeris  Time 
to  other  kinds  of  time  must  be  through  observation  of  these  bodies. 
The  body  easiest  to  observe  is  the  moon,  and  therefore  the  relation 
between  Ephemeris  Time  and  Universal  Time  is  gotten  at  present 
through  the  lunar  motion.  The  principal  techniques  utilised  are 
occultations,  transit  observations,  and  photography  of  the  moon 
against  &  star  background  (the  moon-camera  method  of  Markowitz 
(ref.  32 2J  .  Although  from  the  standpoint  of  operational  definition 
there  is  a  different  kind  of  Ephemeris  Time  for  every  member  of 
the  solar  system,  these  fine  distinctions  are  quite  lost  in  the 
observational  errors  associated  with  measurement  of: 

AT  5  E.T.  -  U.T. 

This  is  true  also  of  the  difference  between  E.T.  as  defined  by 
Brown's  theory  and  BJI.T.  as  defined  by  a  numerically  correct  theory. 
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(7)  The  OT2  time  of  emission  of  the  signal  is  not  the 
same  as  the  time  stated  by  the  signal,  but  may  differ  from  it 
by  many  milliseconds  even  if  the  radio  station  is  operating 
properly.  The  U.  S.  Naval  Observatory  monitors  the  time  signals 
from  a  large  number  of  time-transmitting  stations,  as  do  Greenwich; 
Tokyo  Astronomical  Observatory,  et  al,,  and  the  times  of  receptions 
are  published  as  UT2  times.  Since  the  monitoring  station  whose 
published  reception  times  are  used  in  the  reduction  is  usually 
within  a  hundred  miles  of  the  transmitting  station  whose  signals 
are  used  for  timing  the  occultations,  the  transmission  time  error 
from  this  source  is  certainly  less  then  100/  ...  likewise,  the 
monitoring  errors  themselves  should  be  less  tiian  100// s,  and, 
in  fact,  the  times  of  reception  are  for  most  stations  given  to 
0.1  ms. 


(8)  In  most  cases,  the  occultation  observing  site  is  within 
1*000  km  of  a  time-transmitting  station,  and  one-hop  propagation 
may  be  assumed.  Giving  the  F2  layer  a  height  of  1*00  lqn,  the 
maximum  difference  in  transmission  times  are  given  by: 

Straight-line : 

One-hop : 

Two-hop: 

The  maximum  error  is  therefore  less  than  1  ms,  and  will  in  fact 
usually  be  less  than  0.5  ms. 

(9)  There  are  no  good  theoretical  values  for  the  time 
delay  involved  in  the  movement  of  the  signal  from  the  antenna 
to  the  recorder  input.  The  maximum  delay  occurs  in  the  1000  cps 
UTC  band  pass  filter,  and  direct  measurements  of  the  time  delay 
show  that  this  is  about  1.3  ms. 

7 .  Aberration 


13.3  ms  .  o.3 

13.6  ms _ q . 8 

ll*.l*  ms 


Aberration,  the  apparent  displacement  in  position  of  an 
object,  caused  by  the  finite  velocity  of  light,  is  not  of  any  great 
importance  in  the  equal-limb-line  method.  The  reason  for  this  is 
that  the  observations  are  paired  in  such  a  manner  as  to  cancel  the 
effects  of  aberration  on  the  results.  First,  since  at  a  particular 
point  on  the  earth,  an  occultation  consists  of  the  alinement  of 
star,  lunar  feature,  and  observer,  the  star  and  lunar  aberration 
effects  are  exactly  equal.  Second,  the  time  difference  between 
observations  of  an  occultation  from  two  different  sites  is  so 
small  that  annual  aberration  is  insignificant.  Since  the  interval 
between  the  observations  in  an  equal-limb-line  occultation  is 
almost  always  less  than  2\  hours,  the  difference  in  the  annual 
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aberration  at  the  two  points  in  the  orbit  must  be  leas  than: 

A2  cos  ft  ■  +  k  sin(s-  X)  a  (s-  X) 

A2  (3  *  -  k  sin  ^cos(s-X)  A(s-/l). 

Since  (8-M,  over  a  period  of  3  hours,  is  less  than  3  x  10  the 
aberration  differences  are  less  than  0V06,  and  the  maximum  difference 
is  therefore  less  than  O'.'l.  This  is  insignificant  when  even  a  rough 
approximation  to  the  aberration  correction  is  made;  an  error  of  10  km  in 
estimating  the  amount  of  correction  necessary  would  introduce  less 
than  O'. '00 5  error. 

The  diurnal  aberration,  OV32  cos  (p  is  somewhat  more  important. 

The  changes  in  right  ascension  and  declination  are  given  by 

A  OC  ■  -  0V32  cos  if  cos  H  sec  ^ 

A  S  m  +  0 3 2  cos  (p  sin  H  sin  £ 

(Reference  11).  Latitudes  are  generally  restricted  to  i  60°,  the 
declination  to  ±  30*,  the  hour  angles  H  can  extend  almost  to  ±90*. 

The  diurnal  variation  limits  are  therefore 

Max.  Min. 

A  -  0V37  0 

A  S'-  ±  0V16  0 


From  the  equal-limb-line  view,  the  difference  between  the  values  of 
diurnal  variation  at  the  two  places  is  of  primary  importance.  Here 


Max.  Min. 

A$  .  0  0 

AH  -  85*  0 

A<p  -  ±30*  0 


and  the  variation  in  aberration  is  as  great  as  for  the  single-site  type 
of  observation.  The  error  in  the  estimate  of  the  aberration  is  compounded 
of  the  errors  in 


(1) 

radius  of  the  earth 

(  < 

1:104) 

(2) 

velocity  of  light 

(  < 

1:105) 

(3) 

cosine  of  latitude  of 
observer 

(  < 

1:104). 

It  is  therefore  less  than  OVOOOl  and  can  be  ignored. 
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From  what  is  stated  at  the  beginning  of  this  section,  it 
is  safe  to  neglect  the  errors  arising  from  aberration  effects  in 
the  lunar  orbit.  In  the  Improved  Lunar  Ephemerls,  page  x,  it  is 
mentioned  that  terms  depending  on  the  eccentricity  of  the  lunar 
orbit  are  omitted  (•  periodic  terms  of  0V009  size).  Rules  for 
inclusion  are  given  in  Transactions  of  the  International 
Astronomical  Union.  No,  VII,  page  175  (1950). 


B.  Error  Effects 


At  this  point  it  is  wall  to  summarize  the  conclusions 
of  the  section  just  preceding.  That  section  gave,  as  reasonable 
values  for  the  errors  in  the  various  physical  quantities  affecting 
an  occultation,  the  following: 


Star  Position  Errors 

---- 

<  0'.’8 

Moon  Position  Errors 
Observer's  Position 

Errors 

A0'.'3 

known  datum 

— 

<  l’.'O 

unknown  datum 

-  . 

<60. " 

Time  Error 

— 

<  0?01 

Aberration  Error 

— 

<  0V005 

The  combined  effects  of  these  errors  are  given  by  the  equation: 


4E1 

iE1  9E1  >E15E1  ~b  E1  9E1  ^E1  ^E1 

^  5TT  ***  Til  Ti2 

0  O’-  o. 

AE2 

AE3 

B 

dE2  ^  E2  dE2~dZ2  dE2  dE2  3E2  Tb2 

3Sm  ***  at 

SE3  r>  E3  dE3  9e3  3E3  "5E3  SE3  3E3 

3£m  dE1  ^  E2  d-E3 

111  ,n  *  *  0  O  O 

—  _ 

_ 

A&m 

A1?l 


A%3 

o 

At 


where  the  matrix  elements  have  the  obvious  interpretation.  The 
aberration  error  is  present  in  the  equation,  but  combined  with  the 
star  and  moon  errors. 
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Now  the  quantities  [aE^J  enter  into  the  theory  through' 
the  equation 

[a«tj]  ■  • 

The  are  the  "observed"  differences  in  the  lunar  radii  at 

the  occultation  equal -limb -line  points  for  occultation  j,  while 
the  aj  are  the  corresponding  coefficients.  In  order  that  AE1 
mav  be  determined,  it  is  essential  that  three  or  more  values  of 
be  available.  That  is, 

AE1  =  AEl(  AO"1,  ACT2,  A<T3 . ) 

and  the  equation  is  • 


M  -  [[‘I  Jt[*13}'1 

For  convenience,  let  us  define 


o 


Then  the  successive  derivatives  of  the  (  AE1"),  when  strung  together 
in  a  Taylor's  series,  give 


<f<  a  e1)  -£  .  a.A.-e!.  4„k 

a  uk 


j.k 


^  ^  E1)  3(4^) 

9  (  a  <rJ )  ^  uk 
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and  finally 

[S^E1)] 


This  equation  is  merely  the  mathematical  way  of  saying  that  the 
errors  present  in  the  occultatlon  results  are  to  be  found  from 
the  errors  in  the  radius  of  the  moon  measured  on  the  fundamental 
plane,  magnified  by  inverse  projection  of  these  latter  errors 
from  the  fundamental  plane  back  onto  the  surface  of  the  earth 
(back  into  the  observer's  local  horizon). 


There  are  several  ways  of  continuing  the  solution  of  the 
error-analysis  problem.  One,  which  is  the  simplest  when  there  is 
convenient  access  to  an  automatic  computer,  is  to  compute,  for 
various  combinations  of  values  of  the  u^'s,  values  of  <3~for  unit 
variations  (for  each  combination  of  upvalues)  in  each  of  the 
u^  successively,  i.e.,  compute 

5~ -  6~  (u*,  u2  ,  u3,  u4,  u5,  u®,  u2) 


CT +  ■> 


(TXu  +4u^,  u2,  u 


3  u4 
•  u  » 


u3,  u^,  u2; 


and  so  on,  so  that,  approximately, 
(T+ 

k 

A.  U 


(Because,  as  will  be  remembered  from  the  definitions  in  Part  I, 
is  defined  as 

4<r-  cr  -  <r 

2  1 

where  G~  and  are  the  radii  of  the  moon  as  observed  at 
points  22and  1,  respectively,  the  quantity 

a  A  cr  a<r3  _  a  dr1 

3  “  d  uk  3  uk  * 

The  linearity  of  the  preceding  equations  hence  allows  the  problem 
to  be  worked  for  O'2  and  cr2  separately.  This  separation  is 
assumed  Immediately  above  and  following.) 
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The  Coefficients  hiving  been  evaluated,  it  reaaine  only 
to  coapute,  for  each  k, 


( 


k 


and  to  sua  over  k.  The  result  will,  of  course,  be  applicable 
only  to  the  situation  as  defined  by  the  set  of  u^'s  chosen, 
but  it  would  not  be  too  difficult  te  (a)  repeat  the  procedure 
for  aere  sets  of  u^'s;  (b)  select  a  representative  n*  In  the 
first  place;  or  (c)  select  the  u^  for  aaxiaua  $(AGT)» 

A  representative  case  is  provided  by  an  occultation  of  date 


*1 

1957 

June 

d  h 
14  7h 

25*  0?6774 

119* 

03* 

36V000 

W. 

*1 

34* 

13' 

OOVOOO 

N. 

hl 

0 

aeters 

Star 

28  896 

(AMS  catalog  no.) 

*2 

1957 

June  14^  8^ 

02“  24*000 

*2 

102* 

16 f 

37V679 

W. 

*2 

36" 

12' 

00.585 

N. 

0  : 

aeters 

The  individual  variations  are  given  in  Table  VII.  (The  conversion 
between  and  AX,  A<P  is  gives  closely  enough  by  the  relations 


■  r  cos  ji  A  A 

aean 


-  r 


aean 


A  0 


although  aore  precise  values  can  be  obtained  froa  the  usual 
foraulae  or  froa  one  of  the  well-known  tables  of  functions  on 
the  International  ellipsoid.) 
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Table  VII,  Parameter  Variations 


- 1 

ACri 

A(T2 

4<r, 

-  205.9 

-  205.2 

-  1604.5 

-  1605.4 

212.8 

200.0 

1631.1 

1631.3 

A*i 

-  31.8 

-  31.0 

M 

-  6.4 

0.5 

A-1 

7.3 

5.3 

At 

664.2 

676.0 

1 

*rrm 

265.9 

-  1414.6 

The  rule  for  propagation  of  errors  gives 


where 


(»*d)2 


Atteatloa  oust  be  paid  to  the  aon-equal  time  errors,  but  when  we  are 
through,  we  get,  usimg  the  values  given  at  the  start  of  this  section, 


(s •  d  A  — )  -  23.87295 

partial 

To  this  r.a.s  0"  we  aay  add  the  r.m.s.  error  arising  from  the  aoa- 
equal-liah-liae  observation  circumstances.  If  it  is  assuaed  that 
the  Watts  liab  profiles,  although  good  absolutely  to  only  180  aeters, 
are  good  relatively  to  such  better  -  say  60  meters,  thea  corrections 
can  be  aade  to  oa  this  basis.  That  is,  although  height 


r  .m. 8oC« 

^B, Watts 

"  hB,True> 

100 

meters 

r » in » 8*6* 

^A, Watts 

-  hA  ) 

A,  True'' 

180 

meters 

r.m.s.e 

^B, Watts 

^A,Watts 

30 

me  ter 8 , 

57 

Heacc 


2  -  23.874952  +  3600  meters 

Total 

-  370.0132375+  3600  -  4170.0132375  motor* 

Takiag  aa  bafore  the  worst  case,  we  shall  brlag  this  value  up 
to  the  surface  of  the  earth  for  aa  ebservatloa  at  aa  altitude 
ef  10*  (the  usual  Halt  la  aalactlag  aa  occultatloa  Is  15*). 
Heaca 

d  -  sec2  80*  X  CT 

-  33 

-  33  x  64.578 

-  2131.074.  meters 


be  takea  hare  for  a  short  dlgressloa  to 
"auaerlcal"  error  aaalysle  aad  "algebraic" 
errer  aaalysls.  Uatll  fairly  raceatly,  the  derlvatlea  ef  a 
parameter  ef  a  fuactloa  froa  computed  values  of  the  fuactloa 
at  dlffereat  pelats  la  Its  domala,  rather  thaa  by  computatlea 
froa  aa  aaalytlc  expresaloa  for  the  parameter  has  beea  frewaad 
upoa  by  Mthematlclaas  as  belag  "uaaathematlcal"  aad  less 
desirable.  Thla  attitude,  uafertuaately,  has  beea  adopted  . 
too  eftea  la  the  ether  sdeaces  where  mathematics  Is  applied, 
aad  there  la  very  little  justlflcatioa  say  more  for  the 
attitude.  The  validity  ef  a  method  Is  dstermlmed  largely  by 
Its  cerrectaess  aad  by  Its  ease  of  applicatlea.  Uatll  the  last 
decade,  the  latter  coaslderatlea  made  It  feasible  to  spead  a 
leag  time  derlvlag  a  formula  If  the  formula  could  shortoa  the 
cemputlag  time;  aew,  la  asst  cases,  a  computatlea  caa  bo  carried 
out  much  mere  rapidly  thaa  caa  the  derlvatlea  of  a  formula, 
Realizatloa  Is  gradually  growlag  that  algebra  Is  oaly  a. tool 
to  be  used  la  gettiag  to  aa  ead,  aad  is  aot  aa  oad  la  Itself, 
aad  that  aaalytlc  fermulatieas  must  be  justified  by  the  ead 
for  which  they  are  used,  aet  by  the  cultural  eaviroasmat  la 
which  they  are  used, 

After  this  dlgressloa,  we  caa  go  ea  to  look  at  the 
aaalytlc  foraulatloa  of  tha  errer  aaalysls.  As  was  seea,  the 
major  problem  remalaiag  is  that  of  fiadlag,  for  each  ebservatloa 
point 


A  momaat  will 
discuss  the  roles  of 


58 


This  is  fortunately  quite  simple,  the  basic  equation  to  be 
used  is  that  given  in  Part  1: 

t2  ■  £( l1  -  41  V 

i  v  obs  1  moon' 

a^  -  I  Ha1  -a*1) 

(T  L  ^0  m  0  m 

+  ((2  *  -  A!2)] 

» 0  m  0  m  - 

Die  basic  equations  are 


so  that 


i 

0 

0 

= 

0 

cos 

-  8in  S* 

I3 

0 

sin 

cos 

L  J 

cos  h 
0 

l-s in  h 

1 
0 
0 


0 

1 

0 


0 

ccs  0* 
-  sin  0^ 


sin  h 
0 

cos  h 

0 

sin  0* 
cos  0* 


and  correspondingly 
1 
0 
0 


V 

a 

£l 


cos  f  -  Sin  J 


sin  S, 


cos  f 


1 

0 

0 


0  0 

cos  5  sin  £ 
m 

-sin  cos  $ 


m 


m 


0 

0 
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C08  <  V 

0 

'8in(ATh'  V 


0 

0 

0 

1 

0 


m 


cos  (  Ct  -  0C  ) 
m  * 
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where 


Ik*  Jacebiaa 


la  slaply  derived  frea  these  fonnul&e  since 


if  we  deaete  the  successive  retatlea  Matrices  ia  the  equatieas 
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Expaasieas  ef  particular  tens  are  given  belew: 


3f„ 

IS 

3  t 
To T 

•k 

d<f' 

'•O 

7a 

_ I& 

3  <?' 

_!S 

aA 

3S, 

_2*s 

<5  0, 
a£3 
zv’ 

~&T~ 

3  t 


*  0 


TT  ■  /’<=«<?• 


-f[«*  <?'  sin  £*  +  cea  <?'  ces  S+  ces  ^  ] 
P ces  <3'  sin  sin  h 


sin  k 


-  ^  sin  <f 
P  css  dj1  ces  k 

f  [  ces  Cp  ces  +  sin  sin  5*  ces  k  ] 
P  ces  sin  S*  »i»  k 

/°  [ ces  £*  sin  <?'  -  sin  ces<^  ces  k] 

-  P  ces  S*  ces  q>*  sin  h 

f  [sin  £*  ces  C$r  -  ces  S*  *i®  <? '  ces  k 

-  ^  ces  ces  <£'  sin  h 

-  f3  ces  S*  ces  Q '  sin  h 


ces  h 
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- 0  ■  P  cos  (V '  cos  h 

at  ^  1 


3|o 

“jj“£ —  -  1°  cos  q  •  sia  sia  h 


Similarly,  the  coefficieats  for  the  aooa's  posltioa  errors 
may  be  derived. 


-5*7 

ai 

0 

-  p 

r  ai « 

PiaTT.  L 

a$\. 

- 

- P  \ 

.1.  TT  [ 

sia 

$. 

p  f 

$  . 

sin-rr  Lc” 

< 

P  f 

cos 

£ 

sinTT^  l 

0  a 

*2 

r  r 

sia 

a  a 

a 

sia  rr-B  L 

*  * 

< 

■P  f 

C08 

S. 

3- a. 

sia  TT^  L 

-»<2 

■P  f 

sin 

5. 

<** 

sin  TT^  [ 

arr. 

-  cot  TTm 

*L 

"  C*c  ""« 

tl 

arr. 

-  \  sia  sia  ^  +  cos  cos  ^  cosCc^-  <^) 

a  L 


cos 


sm  Sin(  £Xb  -  <**)] 
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-  ~  f>.~.  I  css  £  cm  $  *!■*( 

•laTT^  L  °*  * 


< .  _jl_  r 

* 

to 

sin  TT",  ]. 

f  r 

.1.  IT.  (. 

*: . 

^L-  1 

SUTT^  | 

a*3 

■  - 

f 

sinir 

■ 

H2 

_ m  • 

ai rm 

-  cot  JTU 

*t3 

■ 

.  f  r 

9  c 

sinrr  L 

■ 

-  ass 

■ 

.  f 

at 

Si  nT, 

+  CM  <T 

>f2. 

.  ^  r 

9  t’ 

sinrr,  l 

£—  f  CM  £  sin  s  -  »i»  £  CM  co,(*«  •  a*>  J 

»  TT\  L  *  m 


•in  £*  cm  S',  -  con  $**  *i»  $n  c..(«  -  <**)j 


£- -  [  CM  S'*  CM  g  8i*(  Q'u  ~  a*)l 

.nTT”  L  M  J 


f  sin  c.s£m-  cm  S*  cm  Sm  coaCoc^  -  <** )j  ^£-a 


,3  d~rr~L 


Sm  CM(  *,  -  «*)  -y^sj  -  «cTT,  “1^ 


la  .  (cm  <5*  CM  £,  +  sin  S*  sin  &„  cm  ( #a  -  %)) 


+  sin 


S,  c.  S.  .!.(«.  -  v  -1^]  -  ~  ^ 
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The  &  and  are  here  assumed  to  be  constant;  as  far  as  the 
three-hour  or  less  interval,  which  is  all  that  concerns  occultations, 
such  an  assumption  is  certainly  valid. 

In  addition  to  the  above  quantities,  there  must  be  taken 
account  of  the  natural  errors  which  creep  directly  into  8 G  through 
our  lack  of  complete  knowledge  of  the  moon's  profile  shape.  Since 
these  errors  are  directly  additive  to  A  (S',  and,  furthermore,  have 
been  discussed  in  preceding  sections,  there  will  be  no  further 
investigation  here  of  the  error,  but  it  will  be  included  after 
evaluation  of  the other  errors.  Note  that  this  same  addition  must 
be  made  to  the  £(AdT)  determined  previously  by  computation. 

Evaluation  of  the  above  expressions  for  typical  values, 
ranges  of  values,  or  for  values  leading  to  a  maximum  for  &{&<?) 
is  not  difficult,  but  is  time  consuming.  The  values  of  the  derivates 
of  the  types 


3  &  ,  d  (k  ,  ^  "n” 

d  t  aTT  y  t 

are  gotten  from  any  lunar  ephemeris.  The  variable  op'  in  these 
equations  is  the  geocentric  latitude,  not  the  geodetic  latitude. 

The  difference  between  the” two  latitudes  is  not  significant  for 
error  evaluation.  Some  simplification  can  be  achieved  by  making 
a  number  of  approximations.  First  of  all,  and  <T#  do  not 
change  appreciably  between  occultations  at  points  1  and  2,  and, 
in  fact,  it  is  customary  in  the  analyBeB  of  the  data  to  use  for 
both  points  an  oL *  and  <5*  computed  for  a  time  half-way  between 
the  times  at  1  and  2.  Hence  we  can  use  the  same  and  for 
1  and  2.  Secondly,  a"*  £_  change  at  rates  of  Ot5  per 
second  and  0316  per  second.  Cver  the  three-hour  maximum  interval 
between  paired  observations,  this  involves  only  about  1?5  change 
lnS^  and  0255  change  in  fa.  Hence  we  can  use  the  same  values  for 
andd^  at  both  ends,  merely  adding  on  a  small  correction  for 
the  change.  Furthermore,  if  a  reasonable  limit  is  being  sought, 
rather  than  an  exact  (but  obviously  not  accurate)  limit,  then  the 
coefficients  of  the  and  $  variations  can  themselves  be  approximated. 
This  is  all  the  more  possible  since  and  $  m  will  not  differ  greatly 
from  and  at  occultation  time.  For  example,  the  correction  to 


is  then 


M 


a  Cl 
a  S, 


-1  [sin  cos  £m  -  cos  y  sin  £m  cos(rtm-<*#)]A  S'* 

siitl^L 

-  COS  y  cos  S'*  sin  (<**  -  <\)  A  <** 

sin-fl^  sin  S'*  cos  <£*  (1  -  cos  (<**  - 
2  r 

-  cos  <4*  sin(  . 

m  m 

6U 


-  -  3/4  X  4.5  X  10-4  ( - 1 1 -  ) 

sin  ~n~~  u 

-  -  3.3  X  10 “4  X  ( - li -  ) 

sin  TTb 


-  -  3.3  X  10“4  X  6.378  X  106  X  57  X  10"6 

—  0.6  meters  per  secern d  of  arc. 

This  quantity  is  thea  te  multiply  the  -A  since  the  ether  parts 
of 


will  cancel  each  other.  The  rest  of  the  differential  corrections 
can  be  gone  through  in  the  same  manner.  The  major  differences  enter 
through  the  portions  which  depend  upon  the  observer’s  coordinates 
since  the  coefficients  here  cannot  be  equated,  the  one  te  the  other. 
The  observer's  coordinates  may  differ  by  30°  to  60°  or  more  in 
latitude;  the  local  hour  angle  h  may  differ  by  90*  or  more  at  the 
two  points.  It  is  to  be  expected  that  the  major  errors  will  creep 
in  at  this  point.  As  can  be  seen  from  the  equations,  the  maximum 
for  any  one  £  t  will  be  about  30  meters  per  second  of  arc. 

The  following  partials ,are  based  on  expected  maximum  variations: 
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PARTIAL  DERIVATIVES  - 


■1 

< 

e 

0 

c 

e* 

ry-7 

< 

S'. 

0 

-  14.67 

17.56 

0 

-  1822.15 

0.08 

Q* 

3.92 

11.10 

-  19.23 

-  1577.94 

-  13.81 

23.92 

<?' 

-  22.21 

21.49 

7.89 

- 

- 

- 

A 

3.92 

11.10 

-  19.23 

- 

- 

- 

t 

3,92 

11.10 

-  19.23 

52.39 

26.87 

-  1.81 

5. 

- 

- 

- 

-  15.95 

1822.24 

0.08 

*■ 

- 

- 

- 

1577.94 

13.81 

-23.92 

7T* 

_ _ 

‘ 

-  1582.16 

-  4.52 

-104390.92 

All  values  expressed  la  aeters/secead  ef  arc. 
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3600» 


where  1^  and  I2  are  the  moments  of  Inertia  of  the  coil  and  pen, 
respectively,  and  and  k^  are  the  torsional  stiffness  of  the 
anchoring  wire  and  pen-coil  wire,  respectively.  The  flux  density 
of  the  (permanent  magnetic  field  is  B,  the  effective  length  of  the 
coil  wire  is  L,  and  the  current  in  the  coil  is  i.  If  the  back 
e.m.f.  of  the  coil  and  drag  are  neglected, 

i  -  E/R 

where  E  is  the  signal  voltage.  To  start  with,  E  will  be  assumed  to 
have  the  form  E-Egl^,  so  that  i  -  i0l^.  The  differential 
equations  of  motion  are  therefore  (references  V~>,  16) 

11  d?  8l _ l^2  (  «2  -  «1>  *  k01  8  '  -Bl1 

dt2  (1) 

12  d2  8;  *  (  Oj  -  «1>  -0 

dt2 

where  O^and  O^are  deflections  of  the  coil  and  of  the  pen,  respective¬ 
ly  from  the  rest,  or  "no  current"  position. 
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Rearranging, 


h  d  A—  -  +  (*oi  ♦  ha5  °1  -  ^2  °2  "  -BLI  (2) 

dtr 

-1=12  »i  *  *2  j!i2_  ♦  kl2  »2  -  0 

dt2 

The  Laplace  transform  of  these  equations  gives, 

*1  p2  01  “  J1  p2  010  +  ^1  +  k12^  01  "  k12  02  “  “BLI  ^ 

-\2  h  *  h.  p2  02  “  h  152  020  +  ki2  02  "  0 

k  "  *01  *  k12 


(Ii  P2  ♦  k)  0!  -  k12  02  -  Ix  P2  0^  =  -BLI 

-k12  0!  ♦  (I2  p2  ♦  k12)  02  -  I2  P2  02O  "  0 

Assume  also  010  -  0  2  020 


Ix  p2  ♦  k  - 

BLI 

‘1 

Then 

02  ' 

_k12 

0 

(5) 

Ix  P2  +  k  - 

k12 

"kl2 

I2  P2 

+  k12 

or 

02  "  • 

-kj.2  BLI 

(6) 

II  I2  P1*  +  P2(I2k  ♦ 

Ilkl2) 

+  (k  -  ki2)k12 

-k12  BU 

m 

TT 

_  (7) 

p1*  +  p2  (L_ 

+  kig) 

CM 

ii 

12 

¥" 
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-k12 


BLI 


II  I2 

+  di  P2  ♦  d2 

The  roots  of  the  expression  in  the  denominator  are 

i  irx  ,  t  ir2, 

so  that  the  transform  of  02  1* 


(8) 


(9) 


.  -k^  BLI 

[1  -  cos  h  rlt 

1  -  cos  h  *2^ 

II  I2 

U  -  ri ' 

1  *i 

r2  / 
r2 

If  the  back  e.m.f .  of  the  coil  is  taken  into  account,  together  with 
the  drag  of  the  pen  across  the  paper,  then  several  new  terms  are 
introduced. 

We  have 


Rj  ♦  BL  d0l 

1  3T 


E 


i  -  l 

R 


BLr  dO^ 
dt 


E 


Eol(t) 


T  *•  .  /  »  ,  rt  -  BLE  B2L2r 

Z1  h  -  k12^°2”°l^  +  ^01®!  *  — h -  - p 


II  «!  -  «l  *  (kol  *  k12>  °1  -  k12®2  • 


BLE 

~TT 


Ij>  ^2  *  n  62  +  k^2  —  0]_) 


(10) 

(11) 

(12) 

(13) 

(14) 
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0  -  .  k12BLE0  [ _ i _ 1  (20) 

R  [(P-ai)(P-a2)(P-a3)(P-aii) 
where  the  aj  are  the  roots  of  the  quartic  In  the  denominator  of 
(19)  hence 
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choice  of  the  values  of  I. .  I_,  k  and  k  removes  the  double 

1*  2*  01  12 

resonance.)  A  somewhat  more  representative  form  for  E  would  have 


been  the  ramp  function. 


71 


APPENDIX  B: 


THE  LIGHT  CURVE 


To  an  observer  looking  at  an  occultation  with  a  telescope,  the  event 
appears  to  be  instantaneous}  in  fact,  it  is  even  possible  for  physiologi¬ 
cal  reasons  for  the  observer  to  "see1*  the  event  before  or  after  it  actually 
occur 8 •  A  recording  of  the  light  intensity,  however,  shows  that  the  event 
is  far  from  instantaneous;  a  noise-free  record  would  show  a  curve  which 
is  similar  to  that  of  a  Fresnel  diffraction  curve  for  a  straight  edge* 

The  slight  (but  important  differences  are  caused  by: 

1)  the  finite  angular  size  of  the  source,  which  makes  the  influence 
of  non-uniform  light  distribution  over  the  star's  surface  felt; 

2)  the  failure  of  the  lunar  surface  to  act  like  a  straight  edge. 
These  two  differences  permit,  in  some  eases,  the  reconstruction  of  the  stel¬ 
lar  structure  or  the  micro-analysis  of  the  form  of  the  lunar  profile 
(References  20,  21).  As  far  a  s  the  geodetic  use  is  concerned,  these  details 
are  of  no  significance  and  are  neglected. 

The  diffraction  pattern  in  its  gross  appearance  can  be  constructed  by 
evaluating  the  Fresnel  integral 


The  variable  v  is  related  to  the  distance  along  the  diffraction  pattern 
by  the  formula 


Evaluating  the  integral  for  A  ■  550  nm,  d  *  380,000  km,  the  intensity  at 
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the  line  observer  at  the  instant  of  the  occultation,  (when  observer, 
noon's  edge,  and  star  are  collinear)  is 

^•obs  "  ^average 

u 

where  Iaverage  t^ie  average  intensity  before  the  immersion  (or 
after  the  emersion).  Prom  the  observer  to  the  first  maximum  is 
approximately  12  meters. 

These  values  must  be  modified  to  take  account  of  the  difference 
between  the  physical  situation  and  the  simplifying  assumptions  which 
were  made  above.  These  differences  include: 

a.  the  superposition  of  diffraction  patterns  from  a  range  of 
wavelengths 

b.  the  different  response  of  the  photomultiplier  tube  at 
different  wavelengths 

c.  the  varying  brightness  of  the  source  at  different  wave¬ 
lengths 

d.  the  filtering  action  of  the  atmosphere  and  optics  on  the 
radiation  from  the  source.  At  a  particular  distance  y  from  the 
diffracting  edge,  the  intensity  as  recorded  would  therefore  be 
represented  by  the  formula 


where  TOO  is  "the  spectral  response  of  the  photomultiplier  tube, 
0(A)  is  the  source  relative  brightness 
H(A)  is  the  atmospheric  filtering  function 
and  K  is  a  scaling  factor. 
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The  RCA  1P21  photomultiplier  tube,  widely  used  in  astronomical 
photometry  and  used  in  the  Amy  Map  Service  occultation  surveys,  has 
an  S-I4  spectral  response*  This  scaled  to  a  value  of  1.00  at  the 
U00  m  maximum,  the  values  are  given  at  50  m  (50  nanometer)  inter¬ 
vals  in  Table  B-l.  The  source  function  can  be  represented  by  the 
Planck  formula  (reference  29) 

h  .  2£l  (e°2/*T  .l)-1 

TV 

where  -  U.992  x  10-1^  erg-cm 

C2  ■  1.U38  cm-deg 

or,  to  a  sufficient  approximation,  by  a  numerical  table  for  the  sun, 
a  0  0  star.  In  Table  B-2  is  given  a  set  of  values  combining  0(A) 
and  H(A),  being  the  solar  radiation  at  sea-level. 

A  simple  picture  of  the  effect  of  these  various  factors  can  be 
gotten  by  looking  back  at  equation  B  (2).  Taking  the  derivative  of 
the  function  y,  and  keeping  the  position  of  the  first  maximum 
constant  at  v  *  1.2  (approximately), 
dy  j  d\ 

y-yp  '  A 

The  effective  range  of  dA  is  about  -  200  around  the  5£0  nm  point, 
so  that  the  maximum  at  different  wavelengths  varies  over  about 
t  2  meters  from  the  550  nm  point.  The  close-by  maximum  and  minimum 
will  be  shifted  by  about  the  same  amount,  so  that  the  diffraction 
pattern  may  be  expected  to  be  considerably  smoothed,  in  the  recording, 
from  what  it  would  be  in  a  monochromatic  recoruing. 


TABLE  B-l 

S-L  RESPONSE 

RESPONSE  («rg/cm2 

A  (na) 

0 

0 

t 

300 

0.0 

350 

0.90 

liOO 

1.00 

L5o 

0.90 

5oo 

0.70 

55o 

O.LO 

600 

0.10 

65o 

0.03 

700 

0.00 

00 

TABLE  B-2 

0 

SOLAR  RADIATION  AT  60°  ZENITH  DISTANCE  AND  AT  SEA  LEVEL 

A  (cm)  x  10^ 

Q(erg  sec“^cnf^)  x  10^ 

QnCqLoo  -  l) 

300 

0 

0 

350 

0.188 

o.Lo 

Loo 

0.L70 

1.006 

L5o 

1.006 

2.1L 

5oo 

1.215 

2.58 

550 

1.190 

2.53 

600 

1.167 

2.L8 

650 

1.173 

2.L9 

700 

1.108 

2.35 

750 

0.867 

1.8L 

800 

0.857 

1.82 
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APPEND H  Cl  ANALYSIS  OF  EQTJAL-LIMB-LINE  OBSERVATIONS 


The  equations  |l3,l$  as  written  enable  one,  apparently,  to  solve 
—  ■) 

for  the  corrections  AEj  at  several  sites  without  the  necessity  of 
making  any  assumptions  about  the  corrections  at  one  of  the  sites.  The 
reason  for  this  apparent  paradox  is  that  the  equations  have  hidden  the 
lunar  orbit  errors  within  the  geodetic  corrections  without  explicitly 
saying  so.  This  is  shown  by  the  following  analyses. 

Equation  (19)  can  be  written  in  vector  form  as  the  difference 
between  two  vectors 


<*j  •  •  (ft-  fj»>2 

M  * 

•  4 

»■  ■  y  —  y  — ^ 

Let  the  true  values  of  /j,  be  Xj,  Xjm  so  that 


Then 

°i 


— V  — ^ 

■  x3  *  **}  ■’ 
f  “  Xjm  +  <^X;jm 


(X2  +  ^X2  "  *2m  "  ^x2m)2 

<r2 


^  — ■? 

Ul  ^AXi  -  X^  -AXijZ 


(X2  -  X2m)2  _  (xl  -  xlm)2  +  (X2  *  X2m) 

<T2  rr~  CT2  ^  ^  ^ 

-  « A.  *  (~X2-X2m)  «  A^2a  +  (Xl-Xlm)  .^Im 

<r2  c~i 


♦  0  (AXj)2 
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The  equation  actually  used  explicitly  includes  only  term  II  without 
specifying  the  disposition  of  terms  I  and  III.  But  term  I  is 
obviously  the  residual  error  in  the  lunar  radius  f*  while 
the  term  III  is  approximately  the  component  of  the  moon's  position 
error  in  the  direction  of  the  tangent  to  the  limb.  Since  term  III 
varies  from  occultation  to  occultation,  it  cannot  be  solved  for, 
and  mu3t  be  absorbed  in  the  overall  error.  Fortunately,  since  the 
two  terms  in  the  pre-factor  are  very  nearly  equal,  differing  in 
direction  only  by  the  libration  in  latitude  and  in  magnitude  by  the 


small  error  in  the  lunar  radius ,  there  is  no  great  error  incurred. 
The  residual  error  is  very  small,  as  follows  from  breaking 


(J-  up  into  components 

J  ?  ■  ^ 

J  jm 


act. 

jx 


where  a<t  is  the  vectorial  error  in  position  of  the  moon's 
jm 
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center, 


is  the  lunar  radius  vector 


A(TX 


is  the  vectorial  error  in  position  of  the  observer. 


Hence 


APPENDIX  D:  SFHEROID-TO-SPHEROID  TRANSFORMATIONS 


A  point  P  with  geodetic  coordinates  (  Xl»  ^1*  **l)  is  given,  the 
coordinates  being  measured  on  and  from  a  spheroid  si  whose  semi -major 
axis  is  ai  and  whose  semi-minor  axis  is  bi.  It  is  required  to  find 
the  geodetic  coordinates  of  P,  (A g*  0g,  hg)  measured  on  and  from 
some  other  spheroid  sg  whose  semi-major  axis  is  ag,  whose  semi-minor 
axis  is  bg,  and  which  is  located  in  some  known  relation  to  si»  The 


position  of  sg  with  respect  to  s^  can  be  defined  by  means  of  two 
rectangular  Cartesian  systems  created  as  follows: 

(1)  The  origin  of  system  tfXij)  is  at  the  center  of  spheroid /s-j) 

fti  M 

(2)  The  I  axis  coincides  with  the  semi-minor  axis  of  fa^ 


and  is  positive  to  the  north. 


(3)  The  f  X^l  axis  lies  in  the  plane  of  the  zero  meridian  of  j si 


and  is  positive  in  the  direction  of  Greenwich.  Then  sg  is  located 
with  respect  to  si  as  follows.  £*lj  i*  translated  parallel  to  itself 


until  the  origins  of  £x^J  and  ^XgJ  coincide.  This  translation 

defines  the  vector  r  ^ 

JA.X^  ,  AXi, 


The  translated 


is  now  rotated  counter-clockwise  about 


until  X^J  lies  in  the  plane  X^  ,  XgJ  ,  defining  the  angle  'f'l. 
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This  new  (translated  and  once-rotated)  system  ^XnJ  is  now  rotated 
counter-clockwise  about  X^  (into  which  X1  went)  until 
coincides  with  X^,  defining  the  angle  4*2*  Finally,  this  translated 
and  twice  rotated  system  is  rotated  counter-clockwise  about  X^  until 
its  remaining  axes  coincide  with  the  corresponding  axes  of  |X2j  , 
defining  an  angle  *3-  Denoting  these  rotations  by  j^RiJ  j^j 
respectively,  we  have 

H  ■ 

H  ■  [RJ{W  -Ml- 

The  solution  to  the  problem  is  then  given  by  the  following  steps  t 

— ^  P(Xi,  xf,  X?) 


(1) 

(2) 


i.  p(Ai>  0i>  hi) 

thru  the  equations 


where 


and 


'A 

(Nq_  +  hx)  cos  0x  cosAi 

A 

■ 

(N1  ♦  thxJcos  0x  sin  A  x 

l3 

1 

|Nx(1  -  «2)thx^  sin  0X 

(3) 


*1 


a^(l  -  e^  sin2  0x)  “2 


®1 


a2  -  b2 
al _ 1 


(U) 


(5) 


2.  P 


3.  P 


[%] 


p  j^j  by  tie  equations 

[«m-  m  <6> 

— >  p(  A2,  02>  h2)  solving  the  equations 
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(7) 


ro*M 
_ ft 

m  m 

(N2  +  h2)  cos  02  cos  A  2 

1 

CVICVi 

1 

■1 

(N2  +  h2)  cos  02  sinA2 
|"n2(1  -  e|)  +  h2|  sin  02 

for  02*  an<*  ^2  ^2  and  e2  ar®  defined  similarly  to  and  e^). 

It  is  step  #3  which  offers  the  greatest  difficulty,  since  it  involves 
the  solution  of  a  system  of  three  simultaneous  non-linear  equations* 
The  longitude  A2  can  b«  found  immediately  from 


^tan  A ; 
I  cot  A . 


\  m  (4/4 
4/4 


(8) 


The  latitude  0g  and  height  h2  may  be  found  by  linearising  the 
equations  to  give  (9) 

o  [-<»,  *  h4)  *  a4^]  sin  cos  A  Q  cos^jCosTy, 

j|3 

♦  hj)  +  JL  cos30jj  sin  0j  sinA0  cos^jsin^j 


1- - 

4 

m 

4 

_  _ 

AA 

A0 

Ah 


and  iterating  to  convergence,  using  A  2,  0i,  h^  for  starting  values* 
Since  the  A,  0,  h  will  be  very  close  together,  the  process  con- 

t 

verges  rapidly.  It  is  sometimes  simpler  to  proceed,  when  an  exact 
solution  is  not  needed,  by  solving  the  equations 
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and  hg  is  confuted  froia  equation  D-7.  While  the  latitudes  0  and 
heights  h  that  result  are  not  exact,  since  0  is  for  an  ellipsoid 
honothetic  to  the  ellipsoid  actually  wanted,  the  difference  in  most 
cases  is  substantially  less  than  1  meter. 

In  practice,  the  angles  will  either  be  identically  zero  by 
definition,  since  it  is  customary  in  determining  "best*  fitting 
ellipsoids  to  enforce  the  parallelism  of  that  ellipsoid's  axes  to 
the  axes  of  the  "best-fitting"  ellipsoid  at  Greenwich,  or  so  small 
that  they  may  be  considered  zero.  Thsipnly  the  translation  vector 
remains.  Furthermore,  in  many  applications,  it  is  merely  a  case  of 
from  one  ellipsoid  to  another,  not  completely  from  one  datum  to 
another.  The  two  ellipsoids  are  then  assumed  tangent  at  the  common 
datum  point,  and  the  transformation  vector  is  given  by 

p  [  ijj  -  p  [iij  -■  jaii]  • 
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APPENDIX  Et  THE  ELLIPSE  OF  ERROR 


A  non-singular  symmetric  matrix  can  be  brought  into  the  diagonal 
canonical  form  by  a  similarity  transformation.  In  the  present  case, 
this  involves  the  diagonalization  of  a  series  of  2  x  2  submatrices 
of  &e3  i  0.  Hence  the  transformation  to  principal  axes  is 


Rl  0  0  ....0 

T 

R!  0 . Ol 

0  R2  . 

0  R2  • 

•  •  • 

•  »  • 

S 

•  • 

•  •  • 

•  •  • 

e  •  • 

•  •  • • »  Rn-1 

•  *  • 

•  •  •  •  •  Rn-1 

—  _ 

L  _ 

where  each  rotation  submatrix  is  of  the  form  (61.2)  or  (61.3),  de¬ 
pending  on  wheter  or  not  the  *  0  or  not. 

The  derivation  of  the  ellipse  of  error  equation  for  the  simple 
case  %  ■  0,  i  ^  2,  is  given  by  Grossmann  (reference  9).  Interpre¬ 
tation  of  the  ellipse  in  this  way  is  not  easy,  however,  but  the 
extension  of  the  ellipse  of  error  to  an  n-dimensional  case  is  quite 
easy.  A  more  general  presentation  is  given  by  Cramer  (reference  10). 

In  dealing  with  the  definition  of  the  n-dimensional  ellipsoid 
of  error,  it  is  simplest  to  consider  the  axes  of  the  ellipsoid  to  be 
defined  by  the  eigen-values  of  the  covariance  matrix,  so  that 


,2 

defines  the  ellipsoid.  here  denotes  the  least-squares 

estimation  of  [aE*"|  ,  and  q2  is  the  parameter  defining  the  size 
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of  the  ellipsoid  (reference  19) «  The  value  of  is  the  same  as  the 
value  of  x^f  this  value  is  chosen  to  represent  the  1-  °(  level  of 
confidence. 

The  confidence  ellipsoid,  in  n-dimensional  space,  is  generally- 
defined  by  the  axes  of  the  coordinate  system  in  vihich  the  measurements 
are  carried  out  and  by  the  standard  deviation  values  gotten  from  the 
variance  matrix.  This  is  not  strictly  analogous  to  the  ellipse  of 
error  used  in  2-dimensional  error  analysis  (reference  y,  Grossman) 

An  extension  of  the  ellipse  of  error  concept  is  desirable,  and  is 
provided  in  the  same  way. 

As  in  the  2-dimensional  case,  a  rotation  matrix,  this  time  of 
3  or  more  dimensions,  is  sought  which  will  convert  the  existing  co¬ 
ordinate  system  to  one  in  which  the  variance  matrix  is  diagonal  and 
for  which  the  variances  are  a  maximum  over  any  other  rotation.  The 
required  rotation  matrix  £  Ry  and  the  corresponding  variance 
matrix  y  are  derived  from  the  initial  variance  matrix  [*■«] 
by  solving  equation 

{[ff-ii]  -  *ii  [Ijj  [Hi]  -0 

for  the  eigen  vectors  QliJ  corresponding  to  the  eigen  values  Z^ 
of  the  original  matrix  |GTyJ  •  The  usual  method  of  computation 
permits  an  easy  ord^irg  of  the  hi  and  rotated  ^  so  that 

0  i  K  j 


li 


l  ii  ti+l  ,  i+1 


^  0 


Hr; 


APPENDIX  F:  TRANSFORMATION  BETWEEN  DATUMS 


The  problem  is  to  convert  corrections  to  the  observation  site  into 
corrections  to  the  corresponding  datum.  Occultations  give  corrections 
|^E*J  to  the  observation  site.  These  corrections  are  made  up  of  two 
portions 

M-[  *]•♦[*][**]« 

The  J  8  arise  from  errors  in  the  actual  survey  connecting  the 
site  to  the  triangulation  system,  while  the  jjd  E^J  arise  from  the 
deflection  of  the  vertical  at  the  datum.  Since  8  andj^E^J  ^ 

enter  into  the  equations  in  different  ways,  the  latter  being  subjected 
to  a  linear  transformation,  it  is  theoretically  possible  to  solve  for 
both.  Practically,  however,  the  number  of  occultations  available  is 
not  sufficient  to  permit  solving  for  a  large  number  of  unknowns,  and 
attention  must  be  focused  on  the|<^EiJ  d,  which  are  the  more  important. 

Ms  can  usually  be  set  equal  to  0,  since  the  occulta tion  sites 
are  connected  to  the  existing  triangulation  by  third-order  survey  and 
the  resulting  error  is  less  than  two  meters  at  the  worst. 

The  formulae  of  Hristow  (reference  8)  permit  adjustment  of  a 
large  number  of  parameters  simultaneously.  They  can  be  written  as 


dA  cos  0 

*11  *12  *13  •  •  •  *16 

dA  cos  0 

d0 

*21  *22  *23  ’  *  *  *26 

d0 

ds 

■ 

*31  *32  *33  *  *  •  *36 

ds 

dA 

•  •  •  •  o  • 

dA 

d<* 

•  •  •  o  •  • 

d  & 

da 

♦  •  o  •  •  • 

da 

. 

Jbs. 

*61  *62  *63  *  *  *  *66 

_ 
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where  s  is  the  geodetic  distance,  A  is  the  azimuth  (from  the  datum), 
£  is  the  flattening,  and  a  is  the  semi -major  axis  of  the  ellipsoid. 
For  the  present  purposes 

dof  »  0 

da  ■  0 


Furthermore,  the  fact  that  P0k6  and  P^  are  connected  through  a  tri- 
angulation  net  implies  that 

ds  ■  0 


While  the  condition  dA  ■  0  is  not  necessary,  it  is  desirable 
because  of  the  paucity  of  occulta tion  data.  Its  enforcement  merely 
keeps  the  internal  orientation  of  the  net  unchanged.  Hence,  the 
transformation  equations  can  be  written  as 


—  — 

A  E1 

■ 

all 

a12 

E1 

A  E2 

—  — 

obs 

a21 

a22 

j 

E2 

datum 


The  individual  elements  of  the  matrix  are 


all 

-  cos  se< 

5 

a12 

=  sec  0S  j  [ 

t(i  2  +  1  ^)aA  + 

(1+t2-  t)2-2t2t)2) 

A0AA  + 

t(l+t2)A02A  -  1 

o 

cos20t(l+t2)AA  -Q  | 

egl 

■  0 

a31 

s  [i-3t(  n 3 

2h2)a^2  - 

-  ^  cos  20  (l+t2+n2)AX2j 


A0  8 


Here 


0oba  “  $  datum 


5  T^obs  “  /^datum 


t  5  tan  0 


e'2  co»2  j2f 
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APPENDIX  G:  THE  ADDITION  OF  CORRECTIONS  FOE  LUNAR  LIMB 
INEQUALITIES 

The  equal-limb-line  method  depends  for  a  great  deal  of  its 
effect  upon  the  assumption  that  at  points  1  and  2  the  occultation 
is  observed,  at  the  same  point  on  the  moon  or  that,  at  any  rate,  in 
the  equation  (the  Taylor  Series  expansion  with  quadrates  and  higher 
terms  dropped) 

Ah  _  J^noon  +  +  _£^oon  A02 

moon  2  c  ^2 

(  P  moon  is  the  lunar  radius  observed  from  point  1;  dt,  a  A  4  0 
'  2  2 

are  the  variations  in  time,  longitude,  and  latitude  causing  the 
occultation  point  on  the  moon  to  change:  and  Ah  is  the 
corresponding  change  in  lunar  height) ,  the  partial  derivatives  are 
small,  This  assumption  may,  in  cases  where  the  quantities  At, 

A  A  2,  A02  ere  large,  but  quite  wrong,  and  a  method  of  correcting 
for  this  error  is  desirable. 

One  such  method  is  to  use  a  chart  of  the  lunar  limb  in  the 
region  being  considered,  and  from  the  chart  to  estimate  the  partials 
in  the  Taylor  Series  linear  terms.  The  definitive  charts  produced 
by  Dr.  C.  B.  Watts  of  the  U.S.  Naval  Observatory  provide  a  test  of 
the  procedure.  These  charts,  with  a  putative  accuracy  of  %  200  meters 
jLn  height,  may  be  used  if  we  assume  that,  although 


h 

1  -  h 

1 

as  200  m 

Watts 

true 

h 

& 

1 

CM 

2 

a 

o 

o 

CM 

Watts 

true 
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yet  VJatts,2  -  '‘..’itts,l  f  htrue,2  -  htrue,l 

consequently,  we  nay  proceed  as  follows* 

u,  For  the  tine  of  occulta tion  at  point-1,  confute  fron  the 
charts  the  lunar  radius  at  the  occultation  point* 

b.  Resolve  the  point-2  shift  in  longitude  and  latitude  into 
components  parallel  to  and  perpendicular  to  the  lunar  profile  at  the 
instant  of  occultation  at  point-1. 

c.  Convert  the  parallel  component  into  limits  which  can  be 
used  on  the  charts,  and  compute  the  height  of  profile  at  the  shifted 
point  on  the  profile. 

d.  Add  the  difference  between  the  two  profiles  to  the  AC  of 
the  survey. 

e.  For  the  predicted  position  angle  corresponding  to  the 
occultation  at  point-2,  and  for  the  profile  at  the  instant  of  occul¬ 
tation  at  point-2  as  seen  from  point-2,  compute  the  lunar  radius* 

f.  Add  the  difference  between  the  two  radii  to  the  ACT 
computed  in  step  d.  above. 

g.  Use  the  new  AO"  to  recompute  the  position  of  site  2. 
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